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This  thesis  deals  with  transport  properties  like  conductance  or  conductivity 
of  systems  in  which  a  small  concentration  of  impurities  has  been  added  to  an 
otherwise  pure  material.  The  small  number  of  impurities  implies  that  the  dis- 
tance between  impurities  is  large  enough  to  neglect  direct  interactions  between 
the  impurities.  Therefore,  all  the  physics  is  a  result  of  the  interaction  of  the 
host  material  with  a  single  impurity.  Depending  on  the  nature  of  the  host  and 
the  nature  of  the  impurity,  the  physics  changes  qualitatively  or  merely  quan- 
titatively. Under  qualitative  changes  I  understand  a  change  in  the  functional, 
often  power-law  dependence  of  physical  properties  on  variables  like  the  temper- 
ature. Quantitative  changes  imply  only  changes  in  certain  constants,  leaving  the 
functional  dependence  unaffected. 

I  will  discuss  two  systems  in  which  qualitative  changes  in  transport  proper- 
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ties  occur.  In  the  first  case,  this  is  due  to  the  special  nature  of  the  impurity.  I 
consider  the  effects  of  a  nonmagnetic  impurity  with  two  states  of  equal  energy  on 
the  conductance  of  clean  metal  point  contacts  and  metal-insulator-metal  tunnel 
junctions.  At  low  temperatures,  the  conductance  shows  a  power-law  dependence 
on  temperature  and  applied  voltage  different  from  the  system  without  the  im- 
purity. Such  behavior  has  been  observed  in  experiments.  I  will  show  that  the 
experimental  and  theoretical  results  are  in  reasonable  agreement. 

The  second  part  of  the  thesis  deals  with  nonmagnetic  impurities  in  a  "d- 
wave"  superconductor,  a  system  probably  closely  related  the  "high-Tc"  super- 
conductors. In  this  case,  it  is  the  "exotic"  host  materials  which  leads  to  effects 
not  anticipated  from  their  "classic"  counterparts.  Again,  the  interplay  of  host 
and  impurity  leads  to  qualitative  changes  in  the  transport  properties.  The  new 
effects  might  help  to  understand  a  puzzle  of  the  low  temperature  properties  of 
the  high-Tc  materials,  for  which  the  standard  model  (which  can  qualitatively 
explain  the  thermodynamic  properties  like  specific  heat)  fails  to  reproduce  the 
experimentally  observed  power-laws  of  the  transport  properties. 
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CHAPTER  1 
INTRODUCTION 

Crystalline  solid  state  materials  can  be  divided  into  two  classes:  pure  materials 
and  materials  with  impurities.  A  pure  material  is  usually  a  lattice  of  one  or  more 
species  of  atoms.  Every  deviation  from  the  perfect  lattice  structure  is  considered 
an  impurity.  Examples  are  interstitional  atoms,  dislocations  and  substitution  of 
some  atoms  by  others  not  present  in  the  perfect  stochiometry.  If  the  concentration 
of  impurities  is  large,  one  approaches  the  amorphous  state,  where  the  long  range 
cristalline  order  is  severely  perturbed.  I  will  not  consider  such  systems. 

I  do  consider  materials  with  small  amounts  of  impurities  such  that  the  dis- 
tance between  impurities  is  large  enough  to  neglect  direct  interactions  between 
impurities.  In  this  case,  the  physics  is  determined  by  the  interaction  of  the  host 
material  with  a  single  impurity.  Depending  on  the  nature  of  the  host  and  the 
impurity,  the  physics  can  change  qualitatively  or  merely  quantitatively.  Often, 
some  physical  properties  are  rather  inert  to  small  amounts  of  impurities,  whereas 
transport  properties  like  conductance  or  conductivity  are  rather  sensitive  and 
show  qualitative  different  behavior  in  the  pure  and  impure  systems. 

In  the  first  part  of  this  thesis  I  consider  the  effects  of  impurities  of  a  spe- 
cial type  on  the  transport  properties  of  clean  metal  point  contacts  and  metal- 
insulator-metal  tunnel  junctions.  At  low  temperatures,  due  to  the  nature  of  the 
impurity,  the  conductance  shows  dramatically  different  behavior  compared  to 
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the  pure  system.  Such  behavior  has  been  observed  in  experiments  by  Ralph  and 
Buhrman  in  1992  [1]. 

The  second  part  of  the  thesis  deals  with  standard,  nonmagnetic  impurities 
in  an  anisotropic,  "d-wave"  superconductor,  probably  close  relatives  to  the  ce- 
ramic, high-Tc  superconductors  discovered  by  Bednorz  and  Miiller  in  1986  [2].  In 
this  case,  it  is  the  "exotic"  host  materials  which  leads  to  effects  not  anticipated 
from  their  "classic"  counterparts,  the  metallic  superconductors  first  observed  by 
Onnes  in  1911  [3].  The  new  effects  might  help  to  understand  the  puzzle  of  low 
temperature  transport  properties. 


CHAPTER  2 

PART  I:  ANDERSON  IMPURITIES  IN  POINT  CONTACTS  AND  TUNNEL 

JUNCTIONS 

Anderson  [4]  introduced  a  model  describing  a  conduction  band  of  electrons 
hybridizing  with  a  localized,  N-fo\d  degenerate  impurity  level  with  on-site  inter- 
action. 

+  5]W^(4cp<,  +  /i.c.)  (2.1) 

P,<7 

This  was  a  step  beyond  the  s-d  exchange  model  introduced  by  Zener  [5]  which 
described  conduction  electrons  interaction  with  a  dilute  concentration  of  localized 
magnetic  moments.  The  Anderson  model  allows  the  study  of  how  such  localized 
moments  are  formed.  It  is  more  physical  and  has  much  richer  behavior,  due  to 
the  interplay  of  the  energy  scales  td,  U  and  T  =  ttNoW^  {No  is  the  density  of 
conduction  electron  states  at  the  Fermi  level).  In  particular,  it  can  be  shown  via 
a  Schrieffer- Wolff  [6]  transformation,  that  in  the  limit  U  =  -2ed  -)•  oo,  VF  oo 
such  that  r/ed  =  NgJ  is  kept  constant,  the  model  maps  onto  the  s-d  model  with 
an  antiferromagnetic  exchange  coupling  J  =  W^/cd. 

H  =  -  I^H<r^pa  -JSYl    (ia^afiCpP  (2.2) 

P>'^  p,p',a,0 
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If  one  lets  U  and  — independently  go  to  infinity,  one  also  creates  a  potential 
scattering  term  due  to  the  lack  of  particle-hole  symmetry.  For  a  wide  conduc- 
tion band  that  is  flat  about  the  Fermi  level,  this  potential  scattering  terms  are 
unimportant  for  the  low  temperature  physics.  However,  they  can  play  important 
roles  for  systems  with  bands  of  vanishing  DOS  at  the  Fermi  level  [7]. 

In  the  following  I  will  show  how  the  formation  of  local  magnetic  moments  and 
their  interaction  with  the  conduction  electrons  as  described  by  the  s-d  model 
can  be  understood  by  looking  at  scattering  processes  of  the  electrons  off  the 
Anderson  impurity.  I  will  then  summarize  the  physics  involved  in  the  screening 
of  these  magnetic  moments  at  low  temperatures,  widely  known  under  the  term 
"Kondo  effect"  [8].  To  obtain  dynamic  properties  I  will  use  a  self-consistent 
approximation  technique  known  as  "Noncrossing  Approximation"  (NCA)  [9-12]. 
The  NCA  can  be  formulated  for  both  equilibrium  and  nonequilibrium  situations 
and  is  especially  suitable  for  the  two-channel  Anderson  model  [13]  in  which  there 
are  two  species  of  conduction  electrons  interacting  with  the  Anderson  impurities. 

I  will  discuss  how  transport  properties  like  the  nonlinear  conductance  can 
be  calculated  for  a  tunnel  junction  with  an  Anderson  impurity  in  the  insulat- 
ing layer  separating  the  metallic  leads.  Similar  expression  can  be  obtained  for 
clean  point  contacts.  Within  the  NCA  I  compute  the  nonlinear  conductance  and 
show  that  there  is  scaling  to  leading  order  for  temperatures  and  applied  bias 
well  below  the  characteristic  energy  scale  of  the  system,  the  Kondo  temperature 
Tk.  This  behavior  is  in  quantitative  agreement  with  experiments  on  clean  metal 
point  contacts  by  Ralph  and  Buhrmann  [1,14],  if  one  assumes  the  presence  of 
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two-channel  Anderson  impurities  in  the  narrow  region  which  defines  the  point 
contact.  Finally,  I  also  discuss  the  effect  of  a  finite  bias  on  the  dynamic  and  static 
susceptibility  for  a  two-channel  Anderson  model.  This  part  of  the  thesis  will  be 
then  summarized  in  the  conclusions. 


CHAPTER  3 

LOCAL  MOMENT  FORMATION  AND  THE  KONDO  EFFECT 


In  this  chapter  I  discuss  how  the  formation  of  local  magnetic  moments  and 
their  interactions  with  the  conduction  electrons  can  be  understood  within  the 
Anderson  model,  described  by  the  hamiltonian  Eq.  2.1.  I  will  then  summarize 
the  screening  of  the  local  moments  at  low  temperatures,  known  as  the  Kondo 
effect.  I  will  also  generalize  to  the  multichannel  Kondo  model  [13]  and  describe 
qualitatitive  differences  of  the  single  and  two  channel  Kondo  models. 

3.1    Formation  of  Local  Magnetic  Moments 

The  Anderson  hamiltonian  Eq.  2.1  describes  two  subsystems  interacting  with 
each  other.  The  noninteracting  conduction  electrons  (described  by  the  operators 
ct,  c)  are  characterized  by  a  wide  band  of  bandwidth  D  and  a  Fermi  energy  ep-  I 
assume  that  the  DOS  N{u)  of  these  electrons  is  flat  around  the  Fermi  levels  for  the 
energies  of  interest  which  are  either  the  temperature  T  or  the  characteristic  energy 
of  the  Kondo  effect,  the  Kondo  temperature  Tk-  For  convenience,  I  measure  all 
energies  from  the  Fermi  level,  that  is  I  take  e/?  =  0. 

The  local  level  has  two  energy  scales:  The  energy  of  a  singly  occupied  level 
and  the  interaction  energy  of  electrons  on  the  local  level  U.  I  will  mainly  consider 
the  case  of  a  two-fold  degenerate  leve,  that  is     =  2.  There  are  4  different  states 
of  the  local  level,  depending  on  its  occupation. 
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•  The  empty  level,  Fig.  3.1  a).  All  energies  are  in  the  conduction  band.  The 
energy  of  this  state  is  zero  (since      =  0). 

•  The  singly  occupied  states,  Fig.  3.1  b).  Since  I  assume  degeneracy,  the 
level  has  the  same  energy  when  occupied  by  either  a  spin  up  or  spin  down 
electron,  namely  e^. 

•  The  doubly  occupied  state.  Fig.  3.1  c).  Due  to  the  Pauli  principle,  the 
only  possible  doubly  occupied  state  is  a  state  with  both  a  spin-up  and  a 
spin  down  on  the  level.  The  energy  of  this  state  is  26^  +  U. 


a) 


b) 


c) 


Figure  3.1:  States  of  an  Anderson  impurity,  a)  The  empty  state  has  no  electron 
on  the  impurity,  all  particles  are  in  the  conduction  band,  indicated  by  the  shaded 
area,  b)  The  singly  occupied  state  has  a  single  electron  on  the  impurity  site. 
In  the  figure,  an  electron  with  spin  down  occupies  the  site.  The  state  with  an 
up-spin  on  the  site  has  identical  energy  (in  zero  magnetic  field),  c)  The  doubly 
occupied  state  has  two  electrons  on  the  impurity  site.  Due  to  the  Pauli  principle, 
the  electrons  must  have  opposite  spin. 
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The  energy  U  is  usually  large  and  positive  because  it  is  an  effective  Coulomb 
repulsion  of  the  electrons  on  the  local  level.  Therefore,  double  occupancy  is 
unlikely,  unless  is  negative  and  je^^l  >  f/.  In  the  following  I  will  assume  \cd\  < 
U/2. 

The  occupancy  of  the  local  level  is  now  determined  by  the  size  and  sign  of  ej. 
I  can  distinguish  three  regimes: 

•  The  empty  impurity  limit,  >>  0.  If  is  large  and  positive,  the  impurity 
level  will  be  essentially  empty,  the  occupancy  is  close  to  zero. 

•  The  mixed  valence  regime,  le^l  ~  0.  Since  ~  ep,  electrons  can  be  in 
the  band  as  well  as  on  the  impurity  without  increasing  their  energy.  The 
occupancy  will  be  about  1/2. 

•  The  Kondo  limit,  \ed\  >  T  =  nNoW^.  If  is  negative  and  but  larger  in 
magnitude  than  the  inverse  lifetime  of  the  singly  occupied  state,  the  level  is 
singly  occupied  essentially  all  the  time  and  the  occupancy  is  close  to  unity. 

It  is  in  the  Kondo  limit  where  the  immediate  interest  lies.  To  see  that  in  this 
case  the  local  level  behaves  like  a  local  magnetic  moment,  consider  the  second 
order  process  of  exchanging  a  spin-up  by  a  spin-down  electron  on  the  level.  Fig. 
3.2. 

In  Fig.  3.2  a)  the  level  is  occupied  by  a  spin-down  electron.  To  empty  the 
level,  one  has  to  use  the  hybridization  matrix  element  W ,  the  electron  has  to  go  to 
the  Fermi  level,  Fig.  3.2  b).  The  local  level  is  now  empty  and  can  be  occupied  by 
either  spin-up  or  spin-down  electron.  The  former  would  restore  the  original  state 
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and  the  process  would  be  ordinary  (nonmagnetic)  scattering.  In  the  latter  case, 
a  spin-up  electron  from  the  Fermi  level  (due  to  energy  conservation)  occupies  the 
local  level,  again  by  use  of  the  hybridization  matrix  element  W.  The  end  result 
is  Fig.  3.2  c),  an  exchange  of  spins  on  the  local  level  (with  an  opposite  change 
in  the  Fermi  sea).  Thus,  this  second  order  process  within  the  Anderson  model 
is  very  similar  to  the  (first  order)  spin  flip  process  of  a  local  magnetic  moment 
of  the  Kondo  model.  The  energy  J  associated  with  this  spin  flip  is  J  =  W'^/ca, 
where  the  denominator  comes  from  the  energy  difference  of  of  the  initial  (flnal) 
state  and  the  (virtual)  intermediate  state.  Observe  that  in  the  Kondo  limit  J 
is  negative,  corresponding  to  antiferromagnetic  coupling  of  the  electrons  to  the 
local  moment. 


iE  Ae  AE 


a)  b)  c) 

Figure  3.2:  Spin-flip  process  in  an  Anderson  model 

In  a  similar  process,  I  could  have  first  let  the  spin-down  electron  hop  on  the 
local  level,  thus  creating  a  doubly  occupied  level,  before  the  spin-up  electron 
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hops  into  the  Fermi  sea.  This  process  creates  a  coupUng  Jy  =  -W^/{t^  +  U). 
In  this  thesis  I  will  only  consider  the  limit  of  large  U,  U  »  \ed\-  This  leads  to 
I  J[;|  <<  I  J|,  so  I  can  neglect  the  processes  involving  doubly  occupied  states. 

The  above  processes  are  only  important  when  the  local  level  is  indeed  occu- 
pied. This  is  the  case  in  the  Kondo  limit,  where  is  large  and  negative.  However, 
if  I  let  \ed\  grow,  the  coupling  J  =  W^/ca  will  become  small.  The  mapping  of  the 
Anderson  model  is  achieved  by  scaling  the  energies  W  and  so  that  J  =  W^/cd 
remains  constant.  In  the  special  case  ea  =  -U/2  oo,  one  obtains  the  Kondo 
model,  i.e.  the  hamiltonian  Eq.  2.2.  If  I  let  first  U  ^  oo  and  then  ->■  -oo,  I 
have  broken  particle-hole  symmetry  already  in  the  Anderson  model,  which  trans- 
lates to  additional  potential  scattering  in  the  corresponding  Kondo  model.  For 
the  flat  density  of  states  A''(a;)  of  the  band  electrons  I  consider,  these  potential 
scattering  terms  do  not  change  the  qualitative  physics  at  low  temperatures. 

3.2    The  Kondo  Effect 

The  Kondo  effect  is  one  of  the  most  fascinating  examples  of  many-body  ef- 
fects, probably  only  topped  by  superconductivity.  There  are  many  books  and 
reviews  about  various  aspects  of  it.  The  fundamental  problem  of  treating  an 
inherently  strongly  interacting  system  has  motivated  many  new  techniques,  of 
which  the  Numerical  Renormalization  Group  (NRG)  developed  by  Wilson  [15,16] 
has  been  the  first  to  give  essentially  exact  statements  [17-22]  about  the  low  tem- 
perature physics  of  the  problem. 

Kondo  [8]  computed  the  self  energy  of  the  Kondo  model  Eq.  2.2  within  third- 
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order  perturbation  theory.  Completely  unexpectedly,  he  found  a  logarithmic  di- 
vergence of  the  third  order  contribution  at  low  temperatures.  The  logarithmic 
contribution  immediately  explained  the  resistivity  minimum  observed  in  certain 
rare  earth  alloys  and  the  logarithmic  upturn  of  the  resistivity  at  temperatures 
below  the  minimum.  However,  the  result  also  implied  that  straightforward  pertur- 
bation theory  would  fail  at  very  low  temperatures,  since  logarithmic  contributions 
were  shown  to  occur  also  at  higher  order  terms,  invalidating  any  theory  which 
would  account  only  for  a  finite  number  of  diagrams.  Abrikosov  [23]  attempted 
to  sum  up  the  leading  order  logarithmic  divergences.  He  found  an  expression  for 
the  self  energy  with  a  divergence  at  a  finite  temperature  ~  T^-.  It  became  ob- 
vious that  nonperturbative  methods  had  to  be  employed  in  order  to  get  definite 
answers  for  the  low  temperature  regime. 

One  approach  was  the  poor  man's  renormalization  group  developed  by  An- 
derson and  Yuval  [24],  which  reinterpreted  the  problems  of  perturbation  theory 
as  a  problem  of  a  increasing  coupling  constant  J{T).  They  showed  that  although 
the  bare  coupling  J  can  be  small  at  high  (room)  temperature,  renormalization 
eflFects  lead  to  a  strong,  at  first  logarithmic  increase  of  J  at  low  temperatures. 
J  increases  until  it  becomes  too  large  for  the  system  to  be  reasonably  described 
by  perturbation  theory.  Finally,  in  a  ground  breaking  work,  Wilson  employed 
his  NRG  to  the  Kondo  model  [15].  He  showed  that  J  indeed  increases  over  any 
bounds  as  the  temperature  is  lowered.  This  implies  a  formation  of  many-body 
bound  state  of  the  impurity  spins  with  the  surrounding  electrons.  Due  to  the 
antiferromagneticity  of  the  coupling,  this  bound  state  is  a  singlet.  For  an  elec- 
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tron  far  away  from  the  impurity,  the  magnetic  moment  is  perfectly  screened. 
Therefore,  it  becomes  conceivable  that  the  low  temperature  state  of  the  system 
is  actually  a  Fermi  liquid,  i.e.  a  system  of  weakly  interacting  fermions.  The  weak 
interaction  comes  from  the  polarization  of  the  screening  cloud  in  the  presence  of 
additional  electrons. 

Depending  on  the  temperature,  the  system  is  in  one  of  three  regimes  which 
are  not  separated  by  phase  transitions  but  have  smooth  crossovers.  In  fact,  the 
intermediate  regime  is  itself  a  crossover  regime  between  the  high  and  the  low 
temperature  regime. 

•  The  high  temperature  regime,  T  »Tk-  The  coupling  J  is  weak.  Electrons 
rarely  scatter  of  the  magnetic  moments.  The  resistivity  is  dominated  by 
electron-phonon  scattering.  The  susceptibility  obeys  a  Curie-Weiss  law 
X  =  7^  with  antiferromagnetic  (positive)  Curie- Weiss  temperature  6. 

•  The  onset  of  screening  at  intermediate  temperatures,  T  ^  Tk-  The  cou- 
pling J  grows  logarithmically,  leading  to  stronger  scattering.  The  screening 
cloud  starts  to  form  but  is  too  incomplete  to  mask  the  growing  coupling. 
Since  the  phonons  are  frozen  out,  the  resistivity  grows  logarithmically.  The 
susceptibility  also  behaves  logarithmically  in  this  regime. 

•  The  low  temperature  regime,  T  «  T^.  The  coupling  J  is  very  strong, 
leading  to  a  quasi-bound  state  of  the  impurity  with  the  surrounding  elec- 
trons. The  impurities  are  screened,  so  that  the  electrons  are  in  a  Fermi  liq- 
uid state.  Consequently,  both  resistivity  and  susceptibility  approach  their 
finite  T  =  0  values  with  quadratic  T-dependence. 
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Nozieres  [25]  presented  a  phenomenological  Fermi  liquid  picture  for  the  low 
temperature  state.  Several  groups  [17-22]  applied  the  Bethe  ansatz  technique 
to  solve  exactly  for  the  energy  spectrum  of  the  Kondo  model.  This  provided 
the  first  exact  analytical  results  for  thermodynamic  properties.  Conformal  Field 
Theory  was  employed  by  Affleck  and  Ludwig  [26,27]  and  Bosonization  techniques 
by  Haldane  [28]  and  Emery  and  Kivelson  [29]  to  obtain  exact  statements  for  self 
energies  in  the  strong  coupling,  low  temperature  regime.  However,  even  today 
there  is  no  technique  available,  which  allows  for  exact  evaluation  of  dynamic 
properties  at  all  temperatures.  The  knowledge  of  the  dynamics  of  the  system  is 
crucial  for  the  determination  of  nonequilibrium  properties.  In  the  next  chapter 
I  will  describe  an  approximative  technique  which  can  provide  qualitative  and 
sometimes  quantitative  correct  results  in  all  regimes. 

3.3    Two  Channel  Kondo  Systems  and  the  Overscreened  Kondo  Effect 

A  basic  assumption  of  the  Kondo  model  Eq.  2.2  is  that  due  to  the  locality 
of  the  interaction  only  the  s-wave  component  of  the  electrons  interacts  with  the 
magnetic  moment.  In  principle,  higher  partial  waves  (p-wave,  d-wave)  could  also 
couple  to  the  local  moment.  If  angular  momentum  is  a  good  quantum  number,  it 
is  conserved.  This  means  that  different  partial  waves  act  like  different  species  of 
electrons  coupled  to  the  local  moments.  This  picture  leads  to  the  multichannel 
Kondo  model  introduced  by  Nozieres  and  Blandin  [13].  In  this  model  M  species 
or  'channels'  of  mutually  noninteracting  electrons  couple  to  a  dilute  concentration 
of  local  moments.  In  general,  the  moments  can  have  any  spin,  but  I  am  mostly 
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interested  in  spin-1/2  impurities.  The  hamiltonian  might  have  the  following 
form: 


H 


(3.1) 


p,aT 


The  new  index  r  =  1 . . .  M  labels  the  M  channels.  *  ' 

Possible  realizations  of  the  model  Eq.  3.1  include  the  quadrupolar  Kondo 
model  introduced  by  Cox  [30],  possibly  relevant  for  certain  heavy  fermion  com- 
pounds [31,32],  and  the  Kondo  model  of  Two  Level  Systems  (TLS)  first  studied 
by  Vladar  and  Zawadowski  [33-37].  In  the  latter  model  there  is  a  confusing 
switch  of  the  active  and  passive  degree  of  freedom.  In  the  straightforward  gener- 
alization of  the  Kondo  model  described  above,  the  impurities  and  the  electrons 
interact  via  a  spin  interaction.  Thus,  the  spin  is  the  active  degree  of  freedom, 
whereas  the  channel  degree  of  freedom  (the  partial  wave  index)  is  unaltered  by 
the  interaction. 

In  the  Kondo  model  of  Two  Level  Systems,  these  roles  are  reversed.  The 
TLS  is  usually  thought  of  as  an  atom  with  two  energetically  degenerate  positions 
embedded  in  a  bulk  matrix  (in  this  case  a  metal).  Since  the  TLS  is  explicitly 
nonlocal,  it  will  certainly  interact  with  different  partial  waves.  On  the  other 
hand,  there  is  nothing  magnetic  in  the  system,  i.e.  one  does  not  have  a  local 
electronic  level  (like  in  the  Anderson  model)  hopping  with  the  atom.  This  means 
spin  is  a  good  quantum  number  and  is  unchanged  by  the  interaction  with  the 
TLS.  However,  the  electrons  can  change  their  partial  waves  index  by  interaction 
with  the  TLS.  This  means  there  are  processes  in  which  the  electron  assists  the 
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TLS  to  change  its  state  (the  atom  hops  from  one  position  to  the  other)  but 
pays  the  price  of  having  to  change  its  partial  wave  index.  Therefore,  the  partial 
wave  index  is  the  active  degree  of  freedom  altered  by  the  interaction  whereas  the 
spin  remains  unchanged  and  serves  as  the  channel  index  in  this  case.  Since  the 
electrons  carry  spin  1/2,  the  above  reasoning  leads  to  a  M  =  2  channel  Kondo 
model,  with  degeneracy  N  =  2  due  to  the  two  states  of  the  TLS. 

The  above  discussion,  although  specific  for  illustrative  purposes,  is  only  an 
example  of  how  to  create  a  channel  Kondo  model  with  a  TLS.  However,  the 
specific  realization  of  the  TLS  and  its  interaction  with  different  partial  wave  or 
superpositions  of  partial  waves  (e.g.  of  different  parity)  is  unimportant  (and 
not  known  for  most  systems,  except  for  certain  glasses).  As  long  as  there  is  an 
interaction  which  allows  changes  in  the  active  degree  of  freedom  (e.g.  partial 
waves,  parity)  and  no  magnetic  interaction,  the  spin  of  the  electrons  provides  the 
two  channels  to  the  Kondo  effect  of  the  active  degree  of  freedom  with  the  TLS. 
This  implies  that  two  channel  Kondo  systems  should  be  rather  common.  The 
question  then  arises  why  systems  with  such  peculiar  behavior  (described  below) 
have  not  been  observed  until  recently. 

A  possible  answer  to  that  question  is  that  the  assumption  of  at  least  near 
degeneracy  of  the  two  states  of  the  TLS  is  not  very  often  fulfilled.  If  the  asym- 
metry A  of  the  two  states  is  too  large,  there  will  be  no  Kondo  effect  at  all  at 
temperatures  T  <  A.  Essentially,  the  system  will  be  weakly  interacting  with 
some  potential  scattering.  At  higher  temperatures  other  eflfects  like  phonons 
will  obscure  any  symptoms  of  two  channel  Kondo  behavior.  Only  if  there  is  a 
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large  enough  range  between  the  temperature  below  which  phonons  are  frozen  out 
and  the  asymmetry  A  can  one  hope  to  observe  the  unusual  two  channel  Kondo 
physics. 

Although  the  hamiltonian  Eq.  3.1  looks  deceptively  similar  to  the  original, 
single  channel  hamiltonian  Eq.  2.2,  the  low  temperature  physics  for  the  case  of 
M  =  2  channels  and  degeneracy  =  2  is  radically  different  from  the  standard, 
single  channel  Kondo  model.  This  can  be  understood  by  noting  that  the  impu- 
rity is  now  overscreened.  In  the  language  of  the  magnetic  Kondo  effect  one  can 
say  that  at  low  T  the  impurity  would  like  to  from  a  singlet  state  with  a  single 
electron.  However,  if  there  are  two  species  of  electrons  (e.g.  red  and  blue  elec- 
trons), each  specie  interacts  with  the  same  strength  with  the  impurity,  so  that  if 
a  red  electron  binds  to  the  impurity,  a  blue  electron  will  do  so  as  well.  Due  to 
the  antiferromagnetic  coupling  the  spins  of  the  red  and  blue  electrons  will  be  of 
opposite  sign  as  that  of  the  impurity.  Because  the  size  of  the  spin  of  the  impurity 
and  each  electron  is  the  same,  the  combined  object  of  impurity  and  the  red  and 
blue  electrons  cannot  form  a  spin  singlet,  but  will  again  be  a  spin  1/2  object  with 
spin  opposite  to  the  original  impurity  spin. 

This  composite  object  will  again  try  to  bind  another  pair  of  red  and  blue 
electrons  with  opposite  spin.  Again,  the  combined  object,  now  consisting  of  the 
impurity  and  two  pairs  of  red  and  blue  electrons  with  opposite  spin,  is  a  spin  1/2 
object  with  the  same  spin  as  the  original  impurity.  This  scenario  can  be  continued 
ad  infinitum.  The  difference  to  the  original  Kondo  effect  is  obvious.  Whereas  in 
the  single  channel  model  the  impurity  spin  was  effectively  screened  at  distances 
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larger  than  ~  vp/Tx  [vp  is  the  Fermi  velocity),  in  the  two  channel  case,  there 
is  no  such  screening.  All  electrons  feel  the  full  strength  of  the  coupling  J.  This 
coupling  does  not  actually  diverge  as  in  the  standard  Kondo  model,  but  it  is  large 
enough  to  radically  change  the  ground  state  of  the  system  at  low  temperatures. 

No  longer  can  the  low  T  physics  be  described  by  a  Fermi  liquid  of  weakly 
interacting  electrons.  All  states  are  many-body  states  of  strongly  interacting 
particles.  Due  to  the  lack  of  screening  the  intermediate  regime  crosses  over  to 
a  new  Non-Fermi  Liquid  (NFL)  regime  at  temperatures  T  «Tk.  Nonpertur- 
bative  methods  are  again  necessary  to  obtain  the  thermodynamic  and  transport 
properties.  All  the  techniques  mentioned  previously  have  been  applied.  The  main 
results  are  as  follows: 

•  A  logarithmic  divergence  of  the  susceptibility  at  low  T,  x  ~  -  log(T/T/<:) 
rather  than  the  finite  Pauli  susceptibility  of  the  single  channel  model. 

•  The  resistivity  shows  nonanalytic  power  law  behavior,  p  =  po{l  -  T'^''^),  in 
contrast  to  the  quadratic  Fermi  liquid  behavior. 

•  The  T  =  0  entropy  is  nonzero,  S  =  1/2  In  2  per  impurity.  This  implies  a 
residual  degree  of  freedom  'half  of  a  spin  1/2  impurity.  The  situation  is 
somewhat  similar  to  a  frustrated  spin  system. 

I  mentioned  above  that  the  Kondo  effect  is  among  the  most  fascinating  exam- 
ples of  many-body  physics.  The  two  channel  Kondo  systems  is  in  some  respects 
even  more  fascinating,  being  a  non-Fermi  liquid  at  low  temperatures.  On  the 
other  hand,  it  is  somewhat  simpler  to  treat  once  one  has  dealt  with  the  strong 
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interactions  precisely  because  there  is  no  'new'  Fermi  liquid  behavior  emerging 
as  one  lowers  T  from  the  intermediate  regime.  It  is  this  feature  which  renders 
the  approximation  I  will  describe  in  the  next  chapter  essentially  exact,  even  and 
in  particular  at  the  lowest  temperatures. 

It  is  rather  amusing  that  the  two  channel  Kondo  model,  which  is  suppos- 
edly more  complicated,  can  be  very  well  described  at  all  relevant  temperatures, 
even  out  of  equilibrium,  by  a  single  technique,  the  Non-Crossing  Approximation 
(NCA).  In  contrast,  for  the  single  channel  Kondo  model,  in  the  focus  of  attention 
for  more  than  30  years,  such  a  technique  has  yet  to  be  found. 


CHAPTER  4 

SLAVE  BOSON  TECHNIQUE  AND  THE  NONCROSSING 
APPROXIMATION  (NCA) 

Before  introducing  the  slave  boson  representation  and  discussing  the  Non 
Crossing  Approximation  (NCA),  I  have  to  say  more  about  the  specific  systems  I 
am  going  to  consider.  The  previous  chapters  dealt  with  dilute  concentrations  of 
Kondo  or  Anderson  impurities  in  a  metallic  bulk.  If  I  apply  a  small  bias  V  to  the 
sample,  the  electrons  about  the  impurities  are  essentially  in  local  equilibrium. 
In  this  case  I  can  apply  linear  response  theory  and  obtain  results  for  e.g.  the 
resistivity. 

I  will  be  mainly  interested  in  nonequilibrium  situations,  where  it  is  no  longer 
possible  to  consider  the  electrons  about  the  impurities  in  local  equilibrium.  Such 
a  situation  can  be  achieved  in  a  tunnel  junction  or  a  point  contact.  For  a  tunnel 
junction  this  quite  obvious:  If  the  impurities  are  located  in  the  insulating  layer 
between  two  clean  metallic  leads,  an  applied  bias  will  drop  almost  entirely  over  the 
range  of  the  insulating  layer,  leading  to  a  nonequilibrium  distribution  of  electrons 
about  the  impurities.  For  a  small  point  contact  the  situation  is  similar.  A  point 
contact  consists  of  two  metal  leads  joined  by  a  narrow  constriction.  Again,  any 
applied  bias  will  drop  over  a  range  of  the  size  of  the  point  contact.  If  this  size 
is  small  and  the  impurities  sit  exclusively  in  the  constriction,  I  again  deal  with  a 
nonequilibrium  distribution  of  electrons  about  the  impurities. 

For  concreteness,  I  will  discuss  the  technique  in  terms  of  a  tunnel  junction. 
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The  model  under  consideration  is  a  multichannel  Kondo  model  as  described  by 
Eq.  3.1,  where  the  impurity  sits  in  the  middle  of  the  insulating  layer  between  the 
two  metallic  leads.  For  the  single  channel  case  such  a  model  was  introduced  by 
Appelbaum  [38]  and  Anderson  [39].  Because  I  assume  the  impurity  in  the  exact 
middle,  the  bias  V  will  only  appear  as  shifts  of  the  chemical  potentials  (i.e.  the 
Fermi  level)  of  the  leads.  This  case  is  generic,  although  more  general  situations 
can  be  considered  and  are  discussed  in  [40]. 

4.1    Slave  Boson  Hamiltonian 

In  chapter  2  and  3  I  have  described  the  equivalence  of  the  Kondo  model  and 
the  Anderson  model  in  the  Kondo  limit.  This  equivalence  can  also  be  used  for  the 
multichannel  models.  Instead  of  trying  to  solve  the  multichannel  Kondo  model, 
I  will  use  the  multichannel  Anderson  model  as  the  starting  point.  The  Anderson 
hamiltonian  of  a  tunnel  junction  with  applied  bias  V  reads 

P,CT,T,Q 

+    E    WMrC^^r  +  h.c),  (4.1) 
p,a,T,a 

where  the  first  term  describes  the  conduction  electron  bands  labeled  by  their 
(pseudo-)  spin  a  =  1 . . .  iV  and  their  channel  index,  r  1 . . .  M.  In  the  presence 
of  an  external  bias  V,  the  conduction  electrons  to  the  left  and  right  of  the  junction 
also  have  different  chemical  potentials  n^,,  a  =  L,R.   The  second  and  third 
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terms  describe  the  impurity  level  and  the  interaction  of  electrons  on  the  level, 
respectively.  The  fourth  term  describes  the  hybridization  of  the  impurity  level 
with  the  conduction  electrons  on  either  side  of  the  junction. 

This  is  a  model  of  an  impurity  level  far  below  the  Fermi  surface  (e^  <<  0) 
hybridizing  with  the  M  channels  of  conduction  electrons.  U  has  to  be  the  largest 
energy  (in  magnitude)  for  the  mapping  of  the  Anderson  to  Kondo  model  (the 
Schrieffer- Wolff  transformation  [6])  to  be  valid.  For  convenience,  I  take  the  limit 
f/  ->■  oo.  This  makes  double  occupancy  of  the  impurity  level  impossible.  Barnes 
[41]  has  introduced  the  'slave  boson  decomposition'  of  the  electron  operators  on 
the  level,  d^^.d^-,,  by  writing 

dl  =  flbr,   d,r  =  blf,  p  (4.2) 

where  /  and  b  are  canonical  fermion  and  boson  operators,  respectively.  The 
constraint  of  no  double  occupancy  in  terms  of  the  electron  operators  d\d  is 
an  inequality,  ZardtrdTr  <  1-  This  means  that  either  the  level  is  empty  or  it  is 
occupied  by  a  single  electron.  With  the  slave  boson  decomposition,  the  constraint 
of  no  double  occupancy  can  now  be  written  as  an  exact  constraint  on  the  new 
operators  /  and  b: 

Q  =  T,fif<r  +  E^blb,  =  i  (4.3) 
In  terms  of  the  slave  boson  operators,  the  hamiltonian  4.1  reads 

p,cr,T,a  T 
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+  E  WMb.o^^r  +  h.c).  (4.4) 

p,(T,T,a 

This  hamiltonian  has  to  be  supplemented  by  the  constraint  Q,  which  restricts 
its  action  to  the  physical  Hibert  space.  The  constraint  can  be  enforced  in  the 
standard  way  by  introducing  a  Lagrange-multiplier  A  and  adding  a  term  \{Q  - 
1).  The  technique  of  how  to  perform  perturbation  theory  and  obtain  physical 
properties  is  described  at  length  in  the  review  of  Bickers  [12].  The  general  method 
involves  computation  in  the  unrestricted  Hilbert  space  where  Q  can  have  any 
value.  Then,  one  projects  onto  the  physical  Hilbert  space  with  Q  =  1  hy  letting 
A  ^  oo  which  eliminates  any  contributions  from  states  with  Q  >  1.  The  states 
with  (5  =  0  can  be  eliminated,  too.  The  result  are  expressions  for  physical 
properties  in  terms  of  the  projected  slave  particle  Green  functions.  Below  I  will 
show  these  expressions  for  the  quantities  of  interest. 

4.2    The  Non  Crossing  Approximation  (NCA) 

4.2.1    Validity  of  the  NCA 

The  strength  of  the  slave  boson  formalism  is  that  it  treats  the  largest  energy, 
the  on-site  repulsion  U  exactly  (for  U  =  oo),  rather  than  perturbatively.  The 
NCA  is  a  self-consistent  perturbation  approximation  for  the  self  energies  of  the 
slave  particles,  E{uj)  (fermion)  and  n(a;)  (boson)  in  the  coupling  of  band  electrons 
to  the  impurity  level  W.  The  second  order  expressions  are  made  self-consistent 
by  inserting  the  dressed  slave  particle  propagators  in  the  Feynman  diagrams 
instead  of  the  bare  propagators  [9-12].    One  can  show  that  this  amounts  to 
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the  summation  of  all  self  energy  diagrams  (in  Matsubara  frequency  space)  for 
which  no  propagators  are  crossing  each  other  (thus  the  name:  Non  Crossing 
Approximation)  [9].  If  one  considers  the  degeneracy  N  as  a  variable  and  performs 
an  l/A'^-expansion,  one  can  also  show  that  the  NCA  accounts  for  all  diagrams 
up  to  order  0{1/N).  All  other  diagrams,  including  all  vertex  corrections,  are 
at  least  of  order  0{1/N^)  and  may  thus  be  neglected  in  the  limit  of  large  N. 
I  am  mostly  concerned  with  degeneracies  N  =  2,  i.e.  Two  Level  systems.  It 
is  not  obvious  that  an  approximation  which  is  valid  for  large  N  gives  at  least 
qualitatively  correct  results  for  N  =  2. 

The  NCA  has  been  very  successful  in  describing  the  single  channel  Kondo 
model,  except  for  the  appearance  of  spurious  nonanalytic  behavior  at  temper- 
atures far  below  the  Kondo  temperature  Tk-  This  does  not  really  limit  the 
application  to  physical  systems,  since  there  is  a  wide  temperature  range  well 
above  the  energy  scale  of  these  NCA-  artifacts  and  still  well  below  the  Tk-  The 
spurious  low-T  properties  are  due  to  the  fact  that  the  NCA  neglects  vertex  correc- 
tions responsible  for  restoring  the  low  T  Fermi  liquid  behavior  of  the  one  channel 
model  [42].  However,  it  has  recently  been  shown  [43]  that  for  the  two  channel 
problem,  where  the  complications  of  the  appearance  of  a  Fermi  liquid  fixed  point 
are  not  present,  the  NCA  does  give  the  exact  low-frequency  power  law  behavior 
of  the  impurity  spectral  function  Ad,{uj)  and  the  suceptibilities  at  zero  T.  There- 
fore, I  expect  to  achieve  a  correct  description  for  quantities  involving  (like  the 
conductance)  and  the  susceptibilities  even  at  the  lowest  temperatures  which  can 
be  reached  within  reasonable  numerical  effort  (about  1/1000  Tk). 
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4.2.2    NCA  for  the  Equilibrium  Case 

The  equations  for  the  self  energies  of  the  retarded  Green  functions  G"'(a;)  = 
(u-Cd-  ^'{co))-^  (fermion)  and  D^{lo)  =  {co  -  U'{uj))-^  (boson)  read  [12,44] 


where  r(w)  =  ■k\W\'^N{u)  =  rN{uj)  {N{lo)  is  the  density  of  states  of  the  band 
electrons)  and  F{uj)  =  1/(1  +  e^'^)  is  the  Fermi  function.  The  self  energies  are 
complex  functions  of  u.  However,  due  to  analyticity  at  finite  temperatures,  the 
real  and  imaginary  parts  are  not  independent.  They  can  be  obtained  from  each 
other  by  means  of  a  Kramers-Kronig  relation,  e.g. 


where  P  in  front  of  the  integral  indicates  Cauchy's  principal  value.  Taking  the 
imaginary  part  of  Eqs.    4.5  and  defining  the  spectral  functions  for  the  slave 


(4.5) 


(4.6) 


particles 


A{uj)  =  -lmG'{uj)/7r 


lmD'{L0)/7T 


I  arrive  at 


A{uj)  M 


J  der{uj  ~  e){l  ~  F{iu  -  e))B{e) 


|G"-(a;)|2  TT 


25 

b(lo)      n  r 

Together  with  Eq.  4.6  (and  the  bosonic  complement)  Eqs.  4.7  form  a  complete  set 
of  equations  to  determine  the  Green  and  spectral  functions  of  the  slave  particles. 

However,  this  is  not  enough  to  compute  the  impurity  spectral  function  Ad{u) 
at  low  temperatures  (the  zero  temperature  limit).  This  becomes  obvious  by 
looking  at  the  corresponding  expression  [11] 

M'-^)  =  ^jdt  e-^'  [  A{t  +  a;)B(e)  +  A{e)B{e  -  lo)]  ,  (4.8) 

where 

Z  =  I  dee-^'[NA{e)  +  MB{e)]  (4.9) 

is  the  partition  function  of  the  impurity  in  the  physical  Hilbert  space.  The 
Boltzmann-factor  does  not  allow  for  a  numerical  evaluation  of  the  integrand 
at  negative  e  if  /?  =  l/kBT  is  large.  It  is  therefore  necessary  to  include  the 
Boltzmann-factor  in  the  spectral  functions  and  find  solutions  for  the  functions 

a{u)  =  e-^^A{u)  ,    b{ij)  =  e'^'^BiLo)   .  (4.10) 


The  corresponding  equations  are  easily  found  from  Eqs.  4.7.  I  can  absorb  the 
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Boltzmann-factor  by  making  use  of  the  relation  F{-uj)  —  1  —  F{uj)  =  e^'^F{u)). 


a(w)  M 


b{uj)  N 


TT 


=  —  /  deT{u)  -  €){F{io  ~  e))b(e) 

TT  J 

I  der{e  -cu){l-  F{e  -  uj)a{e)  .  (4.11) 


With  the  functions  a  and  b  the  equations  for  the  impurity  spectral  function  and 
the  partition  function  are 

Ad{uj)  =  ^  lde[A{e  +  io)b{e)  +  a{e)B{e  -  u)]  (4.12) 


Z  =  j  dt  [Na{e)  +  Mb{e)]  . 


(4.13) 


It  is  instructive  to  realize  that  the  functions  a  and  h  are  proportional  to  the 
Fourier  transform  of  the  lesser  Green  functions,  that  is 


«M  =  ^G<{uj),    G<{t-t')  =  -z{  fit')  fit)) 


=  ^D<iuj),  D<it-t')=i{b\t')bit)) 


(4.14) 
(4.15) 


and  are  proportional  to  the  distribution  functions  of  the  slave  particles.  From 
now  on  I  call  a(a;)  and  6(cj)  the  'lesser'  functions.  Eq.  4.13  tells  me  that  Z  is 
obtained  by  integrating  the  lesser  functions  with  the  corresponding  degeneracies 
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and  M  multiplied.  This  implies  that  Z  is  the  averaged  constraint,  that  is 


(4.16) 


where  the  angular  brackets  denote  the  expectation  value  in  the  general  Hilbert 
space  with  the  projection  on  the  physical  Hilbert  space  imposed  after  the  com- 
putation of  this  expectation  value. 

The  Eqs.  4.6,  4.7  and  4.11  form  a  set  of  equations  which  allow  for  the  construc- 
tion of  the  impurity  spectral  function  (from  which  most  transport  properties 
follow).  The  equations  are  solved  by  iteration  with  gaussians  as  initial  'guesses' 
for  the  spectral  and  lesser  functions  at  high  temperatures.  Once  solutions  for  high 
T  are  found  (that  is,  the  solutions  are  so  stable  under  further  iteration  that  no 
function  value  at  any  frequency  changes  more  than  0.1%),  the  next  job  at  lower  T 
will  use  these  solutions  as  input  functions.  Temperatures  down  to  1/1000  Tk  can 
be  reached  upon  repeating  this  procedure.  After  finding  such  low-T  solutions,  it 
is  favorable  (since  faster)  to  use  these  solutions  as  input. 

A  general  problem  of  iterative  solutions  of  integral  equations  is  convergence 
and/or  convergence  to  the  'right'  (physical)  solutions.  If  I  end  up  with  solutions 
that  do  not  fulfill  the  constraint  Z  =  1,  all  eflfort  was  in  vain.  There  is  an 
elegant  way  to  enforce  the  constraint  at  every  iteration  step  which  I  will  describe 
in  Appendix  A.  Sum  rules  constitute  other  checks  the  solutions  have  to  fulfill. 
The  slave  particle  spectral  functions  should  be  normalized  to  unity,  whereas  the 
impurity  spectral  function  must  fulfill  /fl!e^,(e)  =  1  -  (1  -  l/N)nf  ,  with  nj  the 
fermion  occupation  number.  Also,       =       must  hold  (n^  is  the  real  electron 
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occupation  number  per  spin  (parity)).  All  these  conditions  are  fulfilled  to  within 
0.5%  in  the  worst  cases  (low  T,  etc.)  but  typically  they  do  even  better  (0.1%).  I 
conclude  that  the  solutions  fulfill  all  physical  constraints  very  well. 

A  dramatic  simplification  of  the  above  procedure  can  be  achieved  that  is 
special  to  the  equilibrium  case.  I  note  that  the  relation  between  the  spectral  and 
the  lesser  functions  is  known  and  given  by  Eq.  4.10.  This  motivates  me  to  define 
new  functions  ^(a;)  and  B{uj)  by 


B{co)  = 


F{-u) 


(4.17) 


If  I  could  find  these  functions  by  some  procedure,  the  lesser  and  the  spectral 
functions  would  follow  by  multiplication  of  factors  that  are  well  behaved  in  the 
zero  temperature  limit.  For  the  spectral  functions  the  corresponding  factor  is  by 
definition  F{-uj),  whereas  for  the  lesser  function  I  find 


a{uj)  =  e-^^A{u)  =  e-f^""  F{-uj)A{u)  =  F{oj)A{u) 


(4.18) 


Indeed,  it  poses  no  problem  to  find  equations  for  A{uj)  and  B{lo)  in  very  much 
the  same  way  as  I  found  the  equations  for  the  lesser  functions  from  Eqs.  4.7  for 
the  spectral  functions.  The  results  are 


TT   y  ^  '  F{-UJ) 


|G^(a;)P 


TV 

TT 


(4.19) 
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One  can  convince  oneself  that  the  products  of  Fermi  functions  in  these  equa- 
tions are  completely  well  behaved  in  the  zero  temperature  limit.  Thus,  by  solving 
the  two  Eqs.  4.19  instead  of  the  four  Eqs.  4.7  and  4.11  I  can  save  one  third 
of  the  integrations  I  have  to  perform  per  iteration  (recall  the  two  integrations  in 
Eq.  4.6  and  bosonic  complement).  Using  this  procedure  a  typical  job  on  a  decent 
workstation  takes  between  two  and  five  minutes  (depending  on  T). 

4.2.3    NCA  for  Static  Nonequilibrium 

If  I  apply  a  finite  bias  V  by  setting  /i£  =  +V/2and  hr  =  -V/2,  the  system 
is  no  longer  in  equilibrium.  The  NCA  equations  have  to  be  derived  by  means 
of  standard  nonequilibrium  Green  function  techniques  [45].  Most  important,  one 
can  not  expect  the  simple  relation  Eq.  4.10  (or  some  naive  modification)  between 
the  lesser  and  the  spectral  functions  to  hold.  Therefore,  the  trick  with  introducing 
the  functions  A  and  B  can  not  be  performed.  One  has  to  solve  the  equivalent  of 
Eqs.  4.7  and  4.11  for  the  nonequilibrium  case  without  any  further  simplification. 
This  was  done  first  by  Meir  and  Wingreen  [46].  However,  their  derivation  uses  a 
Lorentzian  density  of  states,  which,  though  formally  favorable,  has  disadvantages 
in  the  numerical  evaluation.  I  follow  an  independent  derivation,  which  allows  for 
an  arbitrary  density  of  states  [40,47]. 

The  NCA  equations  for  static  nonequilibrium  read  (recall  T{uj)  =  TN{uj)) 

A{lo)       M  r 

WW  "^'^^'^  ?  f^"^^"^  -  6  +  /x,)(l  -  F(u;  -  e  +  /x,))] 

B{lo)        N  f  ^ 
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a{uj) 


[  de  b{e)  ^  [T^N{u  -  e  +  fi,){F{uj  -  e  +  /x,))] 


I  de  a{e)  ^  [r,iV(e  -  to  -  fi^){l  -  F{e  -  co  -  ^a)]  (4.21) 


a 


TT 


a 


If  the  density  of  states  A'^(c<;)  were  a  constant,  the  only  difference  between  the 
equilibrium  and  the  nonequilibrium  NCA  equations  would  be  the  replacement  of 
the  Fermi  function  as  a  distribution  function  by  an  effective  distribution  function 


Since  my  density  of  states  is  a  gaussian  with  a  width  larger  than  all  other  energy 
scales  (like  |ed|,  Ftot,  Tk)  this  is  in  fact  the  only  significant  modification  of  the 
NCA  equations  itself.  Numerically,  the  most  crucial  modification  concerns  the 
integration  mesh.  The  proper  choice  of  integration  meshes  is  central  to  the  success 
of  the  iteration  and  is  discussed  extensively  in  Appendix  B. 

The  NCA  can  also  be  generalized  to  the  case  where  both  the  potentials  in 
the  leads  and  on  the  'impurity'  level  are  explicitly  time-dependent,  rather  than 
static.  For  the  case  of  harmonic  oscillations  of  frequency  Q,,  I  have  computed 
the  DC-current  and  DC-conductance  in  the  nonadiabatic  case,  fl  »  Tk  [48]. 
Many  new  features  like  electron  pump  effects  and  side  peaks  in  the  nonlinear 
conductance  are  predicted.  However,  inclusion  of  this  material  does  not  fit  into 
the  frame  of  this  thesis. 


F^Q  given  by  {Ttot  =  F^,  +  F^) 


^eff(^)  =  -  /^l)  +  ^F{e  -  t^^) 

^  tot  L  tot 


(4.22) 


4-3    Current  Formulae,  Conductance  and  Susceptibilities 


31 


4.3.1    Current  Formulae  and  Conductance 

Due  to  the  finite  bias  V  a  current  will  flow  from  left  to  right  through  the 
impurity.  This  current  can  be  computed  from  alone,  if,  and  only  if,  the 
couplings  Ta  of  the  impurity  to  the  leads  are  the  same,  and  if  I  assume  a  wide 
band  that  is  flat  for  energies  of  order  V  about  the  Fermi  level. 

HV)  =  M^^M  I  du;A,{u)[F{u  -  n,)  -  F{u  -  ^r)].  (4.23) 

If  the  couplings  are  not  equal,  another  term  shows  up  [49,50],  involving  the  lesser 
Green  function       of  the  impurity  electrons.       can  be  computed  via 

=  ^  I  dea{e)B{e  -  u)  .  (4.24) 

The  NCA  is  a  current  conserving  approximation  [46].  Therefore,  the  currents 
computed  for  the  left  and  the  right  leads  should  be  the  same  when  formally 
evaluated.  For  general  couplings  and  bands,  the  left  and  right  currents  are  given 
by 

—Ne  f 

h{V)  =  ~^J  duFr^Niu  -  ^il)  -  Aa{u)F{u  -  /i^)]  (4.25) 

Ne  f 

Ir{V)  =  —  j  dwr^Niuj  -  i^n)  [G^{uj)  -  A4u)F{lo  -  ^in)]  (4.26) 

Numerically,  they  agree  better  than  0.5%,  which  sets  a  limit  to  the  uncertainty 
for  the  current  I{V)  =  (7^  +  /r)/2. 

In  order  to  obtain  the  (differential)  conductance  G{V)  =  dI{V)/dV,  I  perform 
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a  numerical  derivative,  using  (/(Vj)  -  I{V2))/{Vi  -  K2),  and  take  it  as  the  value 
of  G  at  the  midpoint  (Vi  +  V2)/2.  This  method  does  not  smooth  fluctuations 
induced  by  the  numerics.  If  the  currents  are  off  by  0.5%,  and  if  the  difference 
[Vi  -  V2)  is  small,  the  corresponding  G  could  fluctuate  wildly  (20%  and  more). 
Such  fluctuations  of  the  conductance  are  not  observed  (at  most  2%  deviations  in 
G),  showing  that  the  currents  have  a  much  smaller  uncertainty  than  indicated 
by  the  difference  of  h  and  Ir  (aside  from  overall  shifts  which  do  not  effect  the 
conductance) . 

To  compute  the  zero  bias  conductance  (ZBC),  I  can  use  either  the  equations 
above  in  the  limit  of  V  ^  0,  or  I  can  use  equations  obtained  from  linear  response 
techniques.  The  latter  method  is  favorable,  since  I  only  have  to  find  the  impurity 
spectral  function  in  equilibrium.  I  use 


C(0,T)  =  iM^|,,^,,(,) 


(4.27) 


for  the  ZBC  in  a  tunnel  junction. 

Within  linear  response,  I  can  also  compute  the  bulk  resistivity  p  of  a  small 
density  of  impurities  in  a  bulk  metal,  p  is  related  to  the  impurity  spectral  function 
via  [12] 

dF{e) 


/uFic) 
de         T{e) , 


(4.28) 


where  r-i(e)  =  CAd{e)  (C  is  impurity  concentration)  is  the  inverse  scattering 
time.  The  constant  in  front  of  the  integral  is  material  dependent.  Observe  that 
because  of  the  "double  inversion"  in  the  above  expression  the  bulk  resistivity  p 
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decreases  with  temperature  above  Tj^.  In  combination  with  the  rise  of  p  at  higher 
temperatures  due  to  phonon  scattering,  this  leads  to  a  resistivity  minimum  at  a 
finite  temperature. 

Most  of  the  calculations  are  done  with  symmetric  couplings,  Tl  =  Tr.  How- 
ever, especially  for  a  tunnel  junction  this  is  probably  not  the  realistic  case.  If 
the  impurity  is  at  some  position  x  within  the  insulating  layer  of  thickness  d,  the 
Tl  will  be  larger  (smaller)  than  r^^,  if  x  is  close  to  zero  (d).  Also,  the  bare 
energy  level  of  the  impurity  will  be  shifted  to  higher  (lower)  values,  if  x  is 
smaller  (larger)  than  d/2.  In  order  to  keep  the  total  coupling  Ttot  =  Tl  +  Tr 
constant  (and  for  simplicity),  I  assume  a  linear  dependence  of  the  T^'s  on  x 
of  the  form  Tl  =  Ttot{l  -  x/d) ,  Fr  =  Ttotx/d.  I  also  modify  according  to 
^diy)  =  ed+{V/2){l-2x/d).  The  latter  modification  turns  out  to  be  insignificant 
as  long  as  V  «  \ed\. 

In  Ref.  [40]  I  show  that  asymmetric  couplings  actually  lead  to  a  conductance 
peak  (zero  bias  anomaly)  that  is  asymmetric  about  zero  bias.  Such  asymmetries 
have  indeed  been  observed  in  experiments  on  tunnel  junctions.  Since  this  is  only 
a  side  issue,  I  refrain  from  just  reproducing  the  data  and  refer  the  reader  to 
Ref.  [40]. 

4.3.2    Tunnel  Junctions  vs.  Point  Contacts 

The  above  formulae  for  the  currents  and  conductances  are  valid  in  a  tunnel 
junction  geometry  where  the  current  must  flow  through  the  impurity.  In  a  point 
contact  two  leads  are  joined  by  a  small  constriction.  A  current  will  flow  through 
the  constriction  without  the  impurity  being  present.  In  fact,  the  impurity  will 
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impede  the  current,  due  to  additional  scattering  in  the  vicinity  of  the  constriction. 
The  question  arises  whether  the  effect  of  the  impurity  in  a  point  contact  is  the 
same  in  magnitude  but  with  a  different  sign  as  in  the  tunnel  junction.  In  Ref.  [40] 
I  consider  a  simple  (one-dimensional)  model.  I  find  that  indeed  the  switch  of  signs 
is  the  only  effect  on  the  current  due  to  the  impurity,  provided  the  transmittance 
t  of  the  constriction  is  close  to  unity,  i  ~  1.  Thus,  in  clean  samples  the  results 
for  the  current  calculated  for  the  tunnel  junction  apply  for  point  contact,  too,  if 
one  subtracts  the  impurity  contribution  from  the  background  current  /„.  If  /„  is 
ohmic,  the  conductance  G{V)  is  shifted  by  a  the  constant  dlg/dV.  Aside  of  this 
shift,  the  conductance  signals  of  a  tunnel  junction  and  a  clean  point  contact  will 
be  the  same  except  for  the  sign. 

4.3.3  Susceptibilities 

The  (dynamic)  susceptibility  is  calculated  using  the  standard  formulae  [11,12] 
from  the  lesser  and  the  spectral  function  of  the  fermions  alone.  The  formula  for 
the  imaginary  part  reads 


The  real  part  can  be  obtained  from  the  imaginary  part  by  means  of  a  Kramers- 
Kronig  relation. 


The  static  susceptibility  Xo  =  Rex{uj  =  0)  follows  from  this  expression. 


(4.29) 


(4.30) 
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In  the  two  channel  Anderson  model  this  susceptibility  is  not  the  magnetic 
susceptibility,  since  the  spin  of  the  electrons  does  not  couple  to  the  'impurity' 
(TLS).  The  magnetic  susceptibility  is  then  given  by  the  formulae  above  with  the 
fermionic  lesser  and  spectral  functions  replaced  by  their  bosonic  counterparts. 


CHAPTER  5 

SCALING  PROPERTIES  OF  SELF  ENERGY  AND  CONDUCTANCE 

In  this  chapter  I  present  the  results  from  the  numerical  evaluation  of  conduc- 
tance and  bulk  resistivity  for  the  case  of  symmetric  couplings  to  the  leads,  using 
the  formulae  discussed  in  the  previous  sections. 

The  results  for  the  two  channel  case  and  the  nonlinear  conductances  have 
been  computed  for  a  point  contact.  For  the  one  channel  case  I  show  the  zero  bias 
conductance  (ZBC)  for  a  tunnel  junction  in  Fig.  5.2  in  order  to  compare  directly 
to  the  bulk  resistivity  p  of  a  metal  with  Kondo-impurities  in  Fig.  5.3.  The  data 
for  the  two  channel  case  have  been  mostly  presented  before  in  Ref.  [47].  The  one 
channel  data  are  presented  to  contrast  the  two  channel  results  and  to  show  the 
failure  of  the  NCA  in  reproducing  the  correct  scaling  exponent  for  A'^  =  2. 

5.1    Linear  Response  Conductances  and  Resistivity 

In  Fig.  5.1  I  show  the  zero  bias  conductance  G{0,T),  for  the  two  channel 
case  {N  =  M  =  2).  As  expected  [26,27],  the  ZBC  shows  T^/^  dependence  at 
low  T  with  deviations  starting  at  about  1/4  Tr-.  Tr-  is  chosen  to  be  the  width 
at  half  maximum  of  the  zero  bias  impurity  spectral  function,  A^,  at  the  lowest 
calculated  T  (see  inset).  The  slope  of  the  T^/^  behavior  defines  the  constant  B^: 

G{0,T)-G{0,0)=B^T'^\  (5.1) 
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Figure  5.1:  Temperature  dependence  of  the  zero  bias  conductance  for  the  M  =  2 
channel  case  in  a  point  contact.  The  zero  bias  conductance  has  T^/^-  dependence 
for  T  <  Tk/4.  This  can  be  used  to  roughly  extract  from  the  experimental 
data.  5s  is  a  material  dependent  constant  which  has  been  divided  out.  Therefore, 
the  slope  of  the  low  T  fit  (dashed  line)  is  equal  unity. 


For  the  one  channel  case  (M  =  1,A^  =  2)  one  would  naively  expect  T^- 
behavior  from  the  exact  solution  of  the  corresponding  Kondo  model  (Fermi  liquid 
behavior  at  low  T).  However,  the  NCA  as  a  large  N  expansion  is  not  able  to 
obtain  this  power  law  for  N  =  2.  Instead,  the  ZBC  shows  dominant  linear  T- 
dependence  at  low  temperatures.  As  discussed  in  the  next  section  the  nonlinear 
conductance  also  has  portions  with  the  corresponding  power  law  (linear  for  A''  - 
2)  as  long  as  T  and  V  are  well  below  T^. 

For  A^  =  4  and  A^  =  6  the  ZBC  has  humps  at  temperatures  below  (Fig. 
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5.2,  for  a  tunnel  junction)  because  the  Kondo  resonance  is  shifted  away  from  the 
Fermi  level  for  N  >  2.  Similar  humps  can  be  seen  in  the  magnetic  susceptibilities 
of  these  systems  [12]. 


Figure  5.2:  Zero  bias  conductance  for  a  tunnel  junction  vs.  temperature  for  the 
M  =  1  channel  model.  The  conductance  for  a  clean  point  contact  would  be 
obtained  by  subtraction  of  this  curve  from  a  (constant)  background  conductance. 
The  graph  for  A'^  =  2  shows  an  almost  linear  T-dependence  at  low  T  whereas  the 
curves  for  spin  degeneracy  N  =  4  and  6  show  nonmonotonic  behavior  (humps) 
The  humps  are  due  to  the  fact  that  the  Kondo  peak  of  the  spectral  function  Ad 
is  shifted  away  from  the  Fermi  energy  by  about  Tk-  For  T  >  Tk  all  the  curves 
fall  like  log{T/Tfc)  for  approximately  one  decade. 

For  a  bulk  Kondo  system  it  is  impossible  to  measure  the  ZBC  of  single  im- 
purities. Instead,  one  can  measure  the  resistivity  p  of  the  bulk  metal.  In  Fig. 
5.3  I  show  the  low  T  parts  of  the  resistivity  for  one  channel  impurities  with 
N  =  2,4, 6.  Only  N  =  6  shows  a  convex  dependence  on  T.  In  fact,  p  behaves  like 
(1  -  aiT/TxY),  as  expected  for  a  Fermi  liquid  [12].  However,  the  NCA  does  not 
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reproduce  Fermi  liquid  behavior  for  iV  =  2.  Again,  this  is  not  surprising,  since 
the  NCA  is  an  expansion  valid  for  large  values  of  A'^. 


Figure  5.3:  Bulk  resistivity  vs.  temperature  for  the  M  =  1  channel  model.  Of 
the  three  curves  only  A^  =  6  has  a  clear  convex  shape  and  falls  roughly  like  at 
low  T.  The  A^  =  2  graph  again  shows  almost  linear  T-dependence.  Note  that  the 
humps  in  the  conductance  for  A^  =  4  and  6  are  not  present  in  the  bulk  resistivity 
P- 


5.2    Nonlinear  Conductance 


Recently,  it  has  been  shown  [47]  that  the  two  channel  model  exhibits  scaling  of 
the  nonlinear  conductance  G{V,  T)  as  a  function  of  bias  V  and  T  of  the  form  [14] 


eV 


G{V,  T)  -  G(0,  T)  =  B^T'^HiiA^))  , 


(5.2) 
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but  there  are  finite  T  corrections  even  for  T  <<  T^r.  Here,  H  is  a  universal  scaling 
function  {H{0)  =  0  and  H{x)  ~  x''  for  a;  >>1)  and  and  A  are  nonuniversal 
constants.  The  exponent  rj  is  1/2  for  the  two  channel  model.  This  scaling  ansatz 
is  motivated  by  the  scaling  of  the  self  energy  of  the  electrons  in  the  variables 
frequency  u  and  temperature  T  as  obtained  by  CFT  in  equilibrium  [26,27].  Since 
such  scaling  is  well  known  to  be  present  in  the  =  2  single  channel  case  [12]  I 
expect  scaling  of  the  conductance  of  the  form  Eq.  5.2,  too,  although  the  exponent 
T]  should  equal  2  (Fermi  liquid  behavior)  in  this  case. 

In  order  to  examine  whether  this  ansatz  is  correct,  the  rescaled  conductance 
is  plotted  as  a  function  of  (eV/kBT)'^.  The  conductance  graphs  for  different  T 
should  collapse  to  a  single  curve  with  a  linear  part  for  not  too  large  arguments, 
e.g.  (eV/kBT)''  <  4  (since  too  large  V  or  T  would  drive  the  system  out  of  the 
scaling  regime).  Such  a  collapse  indeed  happens  for  low  bias  V  <  T.  However, 
for  larger  bias  the  slope  of  the  linear  part  shows  T  dependence  (for  more  details 
see  Ref.  [47]).  This  is  not  contradictory  to  the  scaling  ansatz,  but  it  does  show 
that  there  are  significant  T-dependent  corrections  to  scaling. 

Fig.  5.4  a)  and  b)  show  the  scaling  plots  for  the  cases  M  =  2  and  M  =  1, 
respectively  (A^  =  2  in  both  cases).  Whereas  the  two  channel  case  shows  the 
behavior  described  above  with  the  expected  exponent  rj  =  1/2,  the  NCA  does 
not  give  the  correct  exponent  for  the  one  channel  model,  that  is,  the  standard 
Kondo  model.  In  fact,  the  data  show  scaling,  however,  the  exponent  t]  is  equal 
unity  rather  than  2.  This  seems  to  reflect  the  linear  temperature  dependence  of 
the  conductivity  that  the  NCA  produces  in  this  case.  This  shortcoming  is  another 
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Figure  5.4:  Scaling  plots  of  the  conductance  for  (a)  the  two  channel  case 
(M  =  N  =  2),  and  (b)  the  one  channel  case  (M  =  1,7V  =  2)  for  point  con- 
tacts. With  Fl  =  Tji  and  B-^  determined  from  the  zero  bias  conductance  ,  there 
are  no  adjustable  parameters.  There  are  roughly  two  regimes  in  these  plots. 
For  [eV/kBTf  <  1.5  the  curves  collapse  onto  a  single  curve  and  the  rescaled 
conductance  is  proportional  to  {eV/kBTf.  For  larger  [eV/kBTY  the  rescaled 
conductance  is  linear  on  these  plots.  There  are  substantial  corrections  to  scal- 
ing even  at  T  small  compared  to  Tk-  At  even  larger  biases  this  linear  behavior 
rounds  off,  indicating  the  breakdown  of  scaling.  The  temperatures  are  in  units 
of  Tk,  which  is  different  for  the  two  cases. 
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consequence  of  the  negligence  of  vertex-corrections  within  the  NCA,  or,  in  other 
words,  a  consequence  of  using  an  1/N-expansion  at  A'^  =  2.  On  the  other  hand, 
the  fact  that  scaling  is  present  again  shows  that  the  qualitative  behavior  is  well 
reproduced. 


CHAPTER  6 

SUSCEPTIBILITY  IN  AND  OUT  OF  EQUILIBRIUM 

6.1    Equilibrium  Susceptibility  for  the  Two  Channel  Model 

Finally,  I  also  show  results  for  the  static  and  dynamic  susceptibilities  with  and 
without  finite  bias.  All  data  shown  are  for  the  two  channel  model.  The  results 
for  the  usual  Kondo  model  show  different  power  laws,  but  the  general  behavior 
upon  application  of  a  finite  bias  very  similar. 

In  equilibrium,  in  the  zero  temperature  limit,  the  dynamic  susceptibility  de- 
fined in  Eq.  4.29  is  given  by  a  step  function  of  the  form  [43] 

Imxiuj)  =  cisign(a;)[l  -  C2^Juj/Tk  +  ...  .  (6.1) 

The  NCA  approaches  this  behavior  as  the  temperature  is  reduced.  However, 
as  shown  in  Fig.  6.1  the  step  is  always  broadened  by  the  finite  temperature 
with  the  extrema  located  at  values  which  grow  roughly  with  T^/^.  The  real 
part  follows  by  the  Kramers-Kronig  relation  and  would  diverge  logarithmically, 
but  the  temperature  cuts  off  this  divergence  as  well.  As  a  consequence,  the 
static  susceptibility  Xo  diverges  logarithmically  as  T  approaches  zero,  a  non- 
Fermi  liquid  behavior  predicted  before  [19,26,27,35,43]  and  well  reproduced  by 
the  NCA  technique,  as  shown  in  Fig.  6.2  (circles). 


43 


44 


Figure  6.1:  Dynamic  susceptibility  for  the  M  =  2  channel  case  in  equilibrium 
{V  =  0).  The  graph  shows  how  the  dynamic  susceptibility  approaches  a  leading 
order  step  function  as  the  temperature  goes  to  zero,  as  expected  for  the  two 
channel  case  (see  Eq.  23).  This  is  in  contrast  to  the  linear  behavior  of  the  M  =  1 
channel  model.  The  value  of  the  susceptibility  is  in  arbitrary  units. 

6-2    Nonequilibrium  Susceptibility  for  the  Two  Channel  Model 

Out  of  equilibrium,  the  finite  bias  serves  as  another  low  energy  cutoff,  but 
in  a  nontrivial  manner.  The  extrema  of  the  imaginary  part  of  the  susceptibility 
are  located  at  smaller  absolute  values  than  at  the  corresponding  temperature 
at  zero  bias.  The  logarithmic  divergence  of  the  real  part  is  cut  off  at  about 
V,  so  that  the  static  susceptibility  does  not  diverge  logarithmically  anymore  as 
T  ^  0.  Instead,  it  approaches  a  (K-dependent)  finite  value  with  a  quadratic 
T-dependence.  However,  this  does  not  signal  the  return  of  Fermi  liquid  behavior 
for  T  <V,  since  there  still  is  ,/v/f  behavior  of  the  conductance  for  V  well  below 
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Figure  6.2:  Static  susceptibility  Xo  vs.  temperature  for  zero  and  finite  bias  V. 
In  equilibrium,  Xo  shows  the  characteristic,  expected  logarithmic  divergence  as  T 
approaches  zero  (for  the  two  channel  model).  Out  of  equilibrium,  this  divergence 
is  cut  off  at  a  temperature  somewhat  below  the  bias  V.  The  inset  shows  that  Xo 
falls  with  below  this  cutoff.  For  high  temperatures  T»Tk,  Xo  falls  like  1/T 
(Curie-Weiss  law).  The  value  of  the  susceptibility  is  in  arbitrary  units. 

Tk.  Fig.  6.2  shows  the  T-dependence  of  Xo  for  V  =  I/IOT^  (triangles). 

Similar  behavior  is  observed  in  the  V-dependence.  Now  T  serves  as  a  cutoff 
of  the  logarithmic  divergence.  For  low  V  the  static  susceptibility  saturates  and 
falls  quadratically  with  V.FoxT<V<  Tk  the  susceptibility  then  falls  logarith- 
mically. I  show  the  y-dependence  static  susceptibility  Xo  for  various  T  in  Fig. 
6.3. 

However,  there  is  a  difference  in  T  and  V  in  the  regime  Tk  <  T,V  <  Ftot- 
For  large  T  the  static  susceptibility  behaves  like  1/T,  indicating  Curie-Weiss 
behavior.   However,  for  large  bias  Xo  falls  less  rapidly  Uke  a/V  +  b\og{V)/V. 
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Figure  6.3:  Static  susceptibility  Xo  vs.  bias  V  for  various  temperatures  T.  Xo 
has  a  very  similar  dependence  on  V  and  T  as  long  as  F,  T  <  Tr-  (in  the  scaling 
regime).  Xo  drops  first  like  V"^  and  then  like  \og{V).  However,  for  large  V  »Tk, 
Xo  falls  less  rapidly  with  V  than  with  T,  see  the  next  figure.  The  value  of  the 
susceptibility  is  in  arbitrary  units. 


This  stresses  again  the  different  consequences  of  increasing  T  and  V  once  one  has 
left  the  scaling  regime  T,V  <Tk- 
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Figure  6.4:  Product  of  the  static  susceptibility  Xo  and  temperature  T  (bias  V) 
vs.  T  {V)  on  semi-logarithmic  scale.  The  T-dependence  shows  saturation  at 
high  temperatures  and  therefore  implies  the  Curie  law,  Xo  oc  l/T.  However,  the 
y-dependence  is  linear  at  large  bias,  implying  that  Xo  falls  less  rapidly  with  V 
than  with  T,  Xo  oc  \og{V)/V. 


CHAPTER  7 
CONCLUSIONS  TO  PART  I 


In  conclusion,  I  described  in  detail  the  analytical  foundations  and  numeri- 
cal implementation  of  the  NCA  integral  equations  for  the  one  and  two  channel 
Anderson  model  out  of  equilibrium.  The  algorithms  enabled  me  to  reach  lower 
temperatures  than  previously  obtained,  allowing  the  study  of  new  physics. 

In  linear  response,  I  computed  the  conductance  through  a  tunnel  junction 
as  well  as  the  bulk  resistivity.  The  two  channel  data  for  both  properties  show 
T^/^-behavior  in  agreement  with  results  obtained  by  other  methods.  For  the  one 
channel  model  and  A'^  =  2,  I  find  dominant  linear  behavior  at  low  temperatures. 
For  A^  =  6,  the  bulk  resistivity  drops  with  (Fermi  liquid  behavior),  how- 
ever, the  tunnel  junction  conductance  rises  with  T^,  reaches  a  maximum  below 
the  Kondo  temperature  Tk  and  than  falls  of  logarithmically  at  higher  T.  This 
"hump"  is  associated  with  the  fact  that  the  Kondo  peak  of  the  impurity  spectral 
function  is  shifted  away  from  the  Fermi  level  for  values  of  N  >  2. 

Out  of  equlibrium,  the  nonlinear  conductance  behaves  again  dominant  linearly 
(for  T  <  V  <  Tk)  for  the  one  channel  case  and  N  =  2.  Therefore,  I  can  plot 
the  conductance  as  a  function  of  eV/ksT  and  achieve  scaling  for  modest  bias. 
Whether  similar  scaling  of  the  conductance  but  with  argument  {eV/kBTy  can 
exist  for  the  case  A^  =  6  is  yet  to  be  determined.  The  tunnel  junction  conductance 
falls  with       for  bias  V  <  Tk-  This  is  not  to  reconcile  with  the  hump  in  the 
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T-dependence  of  the  zero  bias  conductance.  If  at  all,  scaling  seems  possible  only 
for  temperatures  well  below  the  temperature  where  the  hump  occurs. 

The  two  channel  data  show  scaling  with  an  argument  {eV/kBTY /2  in  agree- 
ment with  conductance  measurements  on  clean  point  contacts.  It  has  to  be 
pointed  out,  though,  that  this  scaling  as  well  as  the  scaling  in  the  one  channel 
model  is  only  approximate.  Finite  T-corrections  are  observed  in  the  numeri- 
cal data  (but  also  in  the  experimental  data)  for  temperatures  down  to  about 
1/100  Tk. 

I  also  calculate  the  dynamic  and  static  susceptibility  and  discuss  the  modifi- 
cations due  to  a  finite  bias  by  example  of  the  two  channel  model.  The  dynamic 
susceptibility  approaches  a  finite  step  as  T  ^  0,  leading  to  a  logarithmic  diver- 
gence of  the  static  susceptibility  in  this  limit.  A  finite  bias  cuts  off  this  logarithmic 
divergence.  In  a  very  similar  fashion,  the  temperature  cuts  off  the  divergence  as 
the  bias  is  vanishing.  Differences  in  the  bias  and  temperature  dependence  of  the 
static  susceptibility  appear  at  high  bias  and  temperature  (outside  of  the  scaling 
regime). 


CHAPTER  8 

PART  II:  NONMAGNETIC  IMPURITIES  IN  HIGH  T, 
SUPERCONDUCTORS 

For  the  high-Tc  SC's  there  is  mounting  evidence  [51]  that  the  superconducting 
state  involves  pairing  of  electrons  in  a  d-wave  state  rather  than  an  isotropic 
s-wave  state,  Fig.  8.1.  The  most  direct  evidence  is  obtained  through  angle- 
resolved  photo  emission  spectroscopy  (ARPES)  which  shows  a  gap  with  four 
fold  symmetry  and  minima  (nodes  within  the  experimental  resolution)  at  the 
lattice  diagonals  [52].  This  is  in  agreement  with  d-^2_y2  symmetry  of  the  pair 
wave  function. 

Because  this  state  has  nodes  of  the  quasiparticle  excitation  gap  Ajt,  the  density 
of  states  (DOS)  is  in  fact  gapless:  It  grows  linearly  from  the  Fermi  level  and  is 
vanishing  only  right  at  the  Fermi  energy,  see  Fig.  8.2.  This  is  in  contrast  to  a 
classic  SC  with  a  nonvanishing  gap  everywhere  on  the  Fermi  surface,  where  the 
DOS  has  a  gap  2A(T)  about  the  Fermi  energy. 

The  presence  of  nonmagnetic  impurities  in  an  otherwise  pure  d-wave  SC  has 
qualitatively  different  effects  compared  to  their  s-wave  counterparts.  In  classic 
isotropic  s-wave  SC's  nonmagnetic  impurities  (that  is,  pure  potential  scattering) 
leaves  the  SC  completely  unaffected  as  long  as  the  impurity  concentration  is  low 
and  the  scattering  strength  is  small  compared  to  the  Fermi  energy.  This  is  known 
as  Anderson's  theorem  [53]:  instead  of  pairing  between  free  electron  (or  Bloch) 
states  there  is  pairing  of  scattering  states.  Since  the  impurities  are  nonmagnetic 
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a)  s-wave  b)  d-wave 


Figure  8.1:  a)  Isotropic  s-wave,  and  b)  d-c2_j^2-wave  gaps  of  a  superconductor. 
The  solid  lines  show  the  cylindrical  Fermi  surface.  The  dashed  lines  indicate  the 
gaps  in  the  quasiparticle  excitations.  Observe  the  nodes  along  the  diagonals  in 
the  d-wave  state. 

they  cannot  break  the  spin  singlet  pairs  of  such  states.  It  takes  magnetic  impu- 
rities to  achieve  this;  the  spin  flip  scattering  of  electrons  on  magnetic  impurities 
destroys  the  singlet  pairing.  The  consequences  are  e.g.  a  strong  suppression  of 
Tc  upon  doping  with  such  impurities. 

In  a  d-wave  SC  even  nonmagnetic  impurities  are  pair  breaking  due  to  the 
nodes  of  the  gap  on  the  Fermi  surface.  This  has  been  shown  via  self-consistent 
t-matrix  approximation  (SCTMA)  [54, 55]  and  very  recently  also  by  nonpertur- 
bative  methods  [56].  The  effect  on  the  DOS  is  a  further  'filling  in'  of  states  in 
the  pseudogap,  as  shown  by  the  solid  line  in  Fig.  8.2. 

Theories  based  on  this  picture  of  a  d-wave  superconductor  with  nonmagnetic 
impurities  [57,  58]  have  been  qualitatively  successful  [59]  in  explaining  proper- 
ties directly  related  to  the  DOS,  like  the  low  temperature  London  penetration 
depth,  specific  heat  etc.  [60-63]  However,  they  fail  to  describe  the  low  temper- 
ature transport  properties  correctly.  Among  the  more  prominent  examples  for 
this  failure  is  the  low  temperature  behavior  of  the  microwave  conductivity  [64,65]. 
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Figure  8.2:  Density  of  states  (DOS)  of  a  d-wave  superconductor.  The  dashed 
line  shows  the  linear  behavior  of  the  DOS  of  a  pure  superconductor  at  low  ener- 
gies. The  solid  line  shows  the  strong  modification  due  to  nonmagnetic  impurities. 
Observe  that  in  the  'dirty'  case  the  DOS  is  finite  at  the  Fermi  level. 

Experiments  measure  a  linear  T-dependence  for  small  concentrations  of  dopants, 
e.g.  Zn  impurities  in  YBaCuO.  This  is  in  contrast  to  the  quadratic  behavior 
found  in  the  simple  theories  described  above.  It  is  possible  that  this  is  due  to 
the  fact  that  Zn-impurities  in  YBaCuO  are  actually  magnetic,  so  that  the  result 
for  nonmagnetic  impurities  does  not  apply.  However,  the  fact  that  for  a  system 
with  nodes  like  a  d-wave  SC  already  nonmagnetic  impurities  are  pair  breaking 
eliminates  the  big,  qualitative  difference  between  nonmagnetic  and  magnetic  im- 
purities present  in  the  classic  SC's.  One  suspects  that  the  qualitative  difference 
present  in  classical,  fully  gapped  SC's  might  be  only  a  quantitative  one  for  a 
d-wave  SC. 


53 


This  leaves  one  with  a  puzzle,  adding  one  more  to  the  many  involving  the 
high-Tc  materials:  If  the  thermodynamics  is  qualitatively  well  described  by  the 
picture  of  a  'dirty'  d-wave  SC,  how  come  that  the  transport  properties  obtained 
by  the  same  model  and  techniques  fail  to  agree  even  qualitatively  with  exper- 
iments? The  answer  lies  in  the  question  itself:  Obviously  there  are  elements 
lacking  in  the  standard  treatment  when  it  comes  to  describe  two-particle  prop- 
erties involved  in  transport,  whereas  the  single-  particle  properties  determining 
the  thermodynamics  are  properly  described. 

The  objective  of  the  proposed  work  is  to  determine  whether  I  can  achieve  bet- 
ter agreement  with  the  experimental  transport  data  within  a  phenomenological 
model  of  a  d-wave  SC  with  potential  scattering  off  nonmagnetic  impurities  and 
additional  scattering  due  to  self-consistently  determined  order  parameter  pertur- 
bations about  these  impurities.  Such  additional  scattering  can  be  expected  to 
have  small  effect  on  single  particle  properties,  since  the  scattering  is  off-diagonal 
in  particle-hole  space.  Transport  properties,  on  the  other  hand,  are  much  more 
sensitive  to  new  terms  in  the  quasiparticle  scattering  rate  arising  from  this  new 
source  of  scattering.  It  is  therefore  plausible  that  this  mechanism  can  do  the 
trick  of  affecting  only  the  properties  which  are  up  to  now  not  explained  within 
the  simplest  theories  of  a  dirty  d-wave  SC. 

In  the  following  I  will  briefly  review  the  current  evidence  for  d-wave  super- 
conductivity in  the  high-Tg  SC's.  I  will  also  review  the  problems  with  transport 
properties  by  example  of  the  microwave  conductivity.  I  then  introduce  the  stan- 
dard treatment  of  a  pure  and  impure  superconductor,  that  is  BCS-theory  [66] 
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and  the  SCTMA  of  nonmagnetic  impurities.  The  vaUdity  of  this  approximation 
has  been  questioned  [67]  for  a  strictly  two-dimensional  d-wave  SC.  Therefore, 
in  chapter  11  I  discuss  a  nonperturbative  method  to  treat  disorder  in  a  two- 
dimensional  SC  which  shows  that  the  SCTMA  gives  at  least  the  qualitatively 
correct  physics. 

The  main  idea  of  this  part  of  the  thesis  is  the  importance  of  order  parameter 
perturbations  about  the  impurities.  I  will  review  the  results  for  s-wave  SC's 
and  discuss  the  differences  to  order  parameter  perturbations  in  a  d-wave  SC. 
Then  I  will  discuss  the  derivation  of  a  new  T-matrix  for  scattering  off  5-function 
impurities  with  additional  scattering  due  to  order  parameter  perturbations.  The 
knowledge  of  an  analytical  form  of  the  T-matrix  allows  for  the  self-consistent 
determination  of  the  order  parameter  perturbations  SAk{q).  Using  the  SCTMA 
to  deal  with  finite  impurity  concentrations  I  then  obtain  the  density  of  states  and 
the  now  momentum  dependent  scattering  rate.  These  results  will  then  be  used 
to  compute  the  microwave  conductivity.  The  thesis  will  end  with  conclusions. 


CHAPTER  9 

D-WAVE  PHENOMENOLOGY  FOR  THE  HIGH-T,  SUPERCONDUCTORS 


The  evidence  for  pairing  in  a  d-wave  state  (most  likely  dx2_y2)  can  be  roughly 
categorized  in  three  groups 

•  Direct  probes  of  the  order  parameter  itself. 

•  Thermodynamic  properties,  i.e.  properties  directly  derived  from  the  quasi- 
particle  density  of  states  N{lo).. 

•  Transport  properties,  i.e.  two-particle  properties. 

The  first  two  groups  are  reviewed  in  detail  in  Ref.  [51],  the  third  group  is  elabo- 
rately discussed  in  Refs.  [59,64].  I  will  briefly  review  the  main  arguments  of  these 
papers. 

9.1    Direct  Probes  of  the  Order  Parameter  Itself 

ARPES  measures  essentially  the  occupied  part  of  the  angle  and  energy  re- 
solved spectral  function  A{k,tv).  The  measurements  [52]  are  taken  just  above 
Tc  for  the  normal  state  and  well  below  to  measure  the  fully  developed  gap 
in  the  superconducting  state.  The  gap  is  determined  by  the  diflterence  of  the 
"Fermi  level" ,  that  is  the  threshold  between  occupied  and  unoccupied  states,  of 
the  normal  and  superconducting  states.  It  has  turned  out,  that  this  gap  is  angle 
dependent  in  the  bismuth-  (BSSCO)  and  yttrium-  (YBCO)  based  high-T^  SC's, 
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with  maxima  along  the  direction  of  the  lattice  vectors  of  copperoxide  layers  and 
nodes  (as  far  as  the  resolution  of  3  5  meV  of  the  experiments  can  tell)  along 
the  diagonals.  This  is  in  agreement  with  the  gap  of  a  pure  da;2_j,2 -state  and  in 
conflict  with  both  isotropic  or  extended  s-wave  states  (or  combination  of  the 
mentioned  states).  In  principle,  the  finite  resolution  of  the  experiments  allows 
for  an  additional  isotropic  s-wave  component  of  the  order  parameter.  However, 
the  experiments  show  clearly  that  if  such  a  component  exists,  it  must  be  very 
small,  so  that  it  would  have  impact  on  thermodynamic  and  transport  properties 
only  at  extremely  low  temperatures. 

Magnetic  flux  measurements  (for  a  review  see  [51])  are  able  to  determine 
phase  differences  of  the  order  parameter  on  specially  manufactured  Josephson 
junctions.  The  crystals  can  be  grown  such  that  the  grain  boundaries  occur  along 
special  direction  of  the  orthorhombic  lattice,  e.g.  the  a-direction  of  one  grain  has 
a  boundary  with  the  b-direction  of  the  adjacent  grain.  Several  geometries  have 
been  studied  by  several  groups  [68-71].  The  experiments  show  that  there  is  a 
phase  shift  of  tt  between  the  order  parameters  of  the  grains  aligned  as  described 
above.  This  implies  that  there  is  a  relative  minus  sign  of  the  order  parameter 
along  the  a-  and  the  b-directions.  As  a  consequence  magnetic  vortices  have 
been  observed  which  are  quantized  in  halfs  of  the  flux  quanta  (f)o  =  h? /ec.  Such 
vortices  have  actually  been  'seen'  by  imaging  techniques  using  scanning  SQUID 
microscopes.  These  observations  are  again  in  agreement  with  a  di2_y2-pairing 
state,  but  can  not  be  explained  by  neither  an  isotropic  nor  an  extended  s-wave 
order  parameter  (no  relative  minus  sign  between  a-  and  b-  directions). 
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9.2    Thermodynamic  Properties 

The  thermodynamic  properties  are  essentially  determined  by  the  quasi-particle 
density  of  states  N{u)).  As  discussed  in  the  previous  chapter,  the  d-wave  and 
extended  s-wave  pairing  states  show  linear  dependence  for  N{u!)  in  the  pure 
systems,  in  contrast  to  the  fully  gapped  N{uj)  of  an  isotropic  s-wave  state  or  a 
combination  of  a  significant  isotropic  s-wave  component  and  a  d-wave  state.  A 
gap  in  A'^(a;)  wil  lead  to  activated  (exponential)  behavior  of  properties  like  the 
specific  heat  and  the  London  penetration  depth.  However,  gapless  systems  will 
show  characteristic  power  laws. 

Measurements  of  the  London  penetration  depth  [60, 61]  indeed  show  linear 
temperature  dependence  on  clean  samples  as  expected  from  the  pure  d-wave  or 
extended  s-wave  pairing  states,  in  strong  contrast  to  the  exponential  behavior 
of  an  isotropic  s-wave  state.  Even  more  convincing,  upon  doping  with  Zinc  im- 
purities, the  deviation  from  the  T  —  0  London  penetration  depth  of  the  clean 
sample  is  finite  and  increases  quadratically  with  temperature  at  low  tempera- 
tures. This  is  in  complete  agreement  with  the  formation  of  a  finite  N{u  =  0),  as 
expected  from  a  gapless  system  like  d-wave  or  extended  s-wave,  in  the  presence 
of  nonmagnetic  impurities  (see  chapter  above).  An  isotropic  s-wave  SC  would 
be  unaffected  my  nonmagnetic  impurities  (Anderson's  theorem),  so  a  change  of 
power  law  is  not  comprehensible  with  this  pairing  state  (even  if  one  somehow 
could  explain  the  power  law  for  the  pure  system  in  the  first  hand). 

Specific  heat  measurements  have  been  more  difficult  because  of  the  large 
phononic  contribution  that  dwarfs  the  electronic  part  everywhere.  Only  at  the 
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phase  transition  does  the  electronic  show  up  as  a  minor  blimp  added  to  the  over- 
all increase  due  to  the  phonons.  Recently,  several  experiments  [62, 63]  with 
twinned  and  untwinned  crystals  have  been  able  to  subtract  the  phononic  contri- 
bution by  careful  comparison  to  reference  samples.  They  claim  that  there  are 
both  linear  and  quadratic  contributions  to  the  specific  heat  C  in  zero  magnetic 
field.  In  an  applied  field  H,  C  (x  H^l"^.  Both  the  quadratic  temperature  depen- 
dence and  the  square  root  dependence  on  the  magnetic  field  are  in  nice  agreement 
with  an  order  parameter  with  nodes,  like  a  d-wave  or  extended  s-wave  SC  [72]. 
An  isotropic  s-wave  SC,  in  the  other  hand,  will  show  activated  (exponentially 
suppressed)  behavior  in  the  T-dependence,  due  to  the  finite  gap  in  the  quasipar- 
ticle  DOS.  For  such  a  SC  the  magnetic  contribution  to  the  specific  heat  comes 
excitations  localized  in  the  core  of  magnetic  vortices.  The  number  of  the  vortices 
is  oc  H,  therefore  the  resulting  specific  heat  would  be  also  oc  H,  in  contrast  to 
the  experiments. 

The  linear  frequency  dependence  of  the  DOS  of  the  order  parameters  with 
nodes  can  explain  the  quadratic  contribution  but  fails  to  explain  the  linear  con- 
tribution in  the  very  clean  (nominally  pure)  samples.  The  fact  that  untwinning 
of  the  samples  strongly  eflFects  the  coefficient  of  the  linear  term  but  barely  effects 
the  quadratic  term  or  the  field  dependence  suggests  that  the  linear  term  might 
not  be  related  to  the  question  of  order  parameter  symmetry  at  all.  There  are 
speculations  that  the  linear  term  could  be  due  to  Two-Level-Systems  located  be- 
tween the  copperoxide  planes,  which  would  probably  not  affect  the  experiments 
probing  the  in-plane  penetration  depth  or  transport.  Still,  the  linear  term  is  an 
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open  question  that  needs  to  be  resolved  before  the  case  can  be  closed  in  favor  of 
the  order  parameters  with  nodes. 

I  conclude  that  the  thermodynamics  of  pure  and  impure  high-Tg  SC's  follows 
naturally  from  the  assumption  of  a  d-wave  or  extended  s-wave  pairing  state, 
but  is  in  stark  contrast  to  isotropic  s-wave  pairing.  However,  thermodynamic 
properties  cannot  distinguish  between  d-wave  or  extended  s-wave  pairing  due  to 
the  identical  density  of  states  N{uj). 

9.3    Transport  Properties 

Transport  properties  are  the  most  indirect  measure  of  the  order  parameter. 
Nevertheless,  because  of  fact  that  they  also  probe  the  relaxation  time  in  addition 
to  the  DOS  they  also  have  the  potential  to  bring  out  the  differences  between 
otherwise  similar  states.  As  with  the  thermodynamic  properties,  studies  [59, 
64, 65]  of  the  microwave  conductivity  a{T)  of  SCs  with  nonmagnetic  impurities 
have  been  able  to  obtain  rough  agreement  with  an  order  parameter  with  nodes 
of  the  d-wave  or  extended  s-wave  type.  These  studies  also  have  shown  that 
admixtures  of  an  isotropic  s-wave  component  must  be  very  small  to  be  permitted 
by  the  experimental  data.  However,  especially  the  low  temperature  data  are 
not  in  good  qualitative  agreement  even  with  a  d-wave  SC.  a{T)  rises  linearly 
with  temperature  in  experiments  with  small  doping  of  Zc-impurities,  rather  than 
quadratic  as  the  simplest  theories  for  a{T)  in  a  dirty  d-wave  SC  predict  [59]. 
Although  these  theories  use  a  phenomenological  model  of  inelastic  scattering  to 
account  for  the  downturn  of  a{T)  at  temperatures  closer  to  Tc,  it  is  hard  to 
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believe  that  the  inaccuracy  in  the  treatment  of  inelastic  scattering  could  account 
for  discrepancies  in  the  low  temperature  regime. 

It  is  therefore  reasonable  to  assume  that  the  standard  treatment  of  the  elastic 
scattering  off  the  nonmagnetic  impurities  is  not  quite  right.  I  have  investigated 
the  possibility  that  additional  scattering  due  to  impurity-induced  order  parame- 
ter perturbations  can  have  significant  impact  for  the  transport  properties  in  SC's 
with  nodes. 


CHAPTER  10 

BCS-HAMILTONIAN  AND  T-MATRIX  FORMULATION 


The  theory  of  impurities  in  a  superconductor  rests  on  two  cornerstones  of 
condensed  matter  theory,  namely  the  BCS-theory  [66]  of  superconductivity  and 
the  self-consistent  Green  function  treatment  of  impurity  scattering  introduced 
by  Abrikosov  and  Gor'kov  [73].  To  clarify  notation  and  as  an  introduction  I  will 
give  a  brief  overview  of  the  physics  involved. 

10.1    BCS-Theory  of  Superconductivity 

In  the  classic  SC's  the  formation  of  pairs  of  electrons  (the  Cooper  pairs)  is  a 
result  of  an  instability  of  the  Fermi  sea  in  the  presence  of  the  retarded  electron- 
phonon  interaction.  If  the  interaction  is  weak,  BCS-theory  (a  mean  field  theory) 
gives  a  good  description  of  the  physics.  In  high-r^  SC's,  although  there  are 
plenty  of  suggestions,  the  underlying  pairing  mechanism  is  not  known.  However, 
consensus  has  it  that  phonons  play  only  a  minor  role  in  bringing  about  super- 
conductivity. Spin  fluctuations  but  also  mechanisms  involving  charge  degrees  of 
freedom  are  seriously  considered.  In  these  cases  the  basic  energy  scale  ujo  is  of 
the  order  of  the  Fermi  energy  e^,  rather  than  the  Debye  frequency  ujd  involved 
in  phononic  mechanisms.  This  change  in  energy  scale  allows  (in  principle)  for  the 
high  Tc's  observed  without  the  necessity  of  strong  pair  interactions  V,  since  within 
BCS  Tc  =  \.Uujoex^{-l/\NoV\).  (It  is  questionable,  though,  that  ujo  and  V  are 
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independent  parameters  in  reality,  as  experience  with  phonon  mediated  supercon- 
ductivity shows.)  It  is  therefore  plausible  to  use  the  weak-coupling  BCS-theory 
rather  than  the  much  more  involved  strong  coupling  Eliashberg-theory  [74]  for 
the  description  of  superconductivity,  at  least  as  a  starting  point. 

BCS-theory  is  based  on  the  BCS  hamiltonian,  a  hamiltonian  in  which  a  four- 
fermion  interaction  has  been  replaced  by  the  interaction  of  the  quasiparticles 
with  a  self-consistently  determined  mean  field  A^,  the  order  parameter.  It  is 
convenient  to  use  the  Nambu-spinor  [75]  =  (cl  .^,c^k,i)  and  its  hermitian 
conjugate  in  order  to  simplify  the  treatment  of  both  (quasi)-particle  and  hole 
degrees  of  freedom.  The  c-  operators  are  standard  second  quantized  fermionic 
operators.  In  terms  of  the  Nambu-spinors  the  BCS-hamiltonian  of  a  pure,  bulk 
superconductor  reads  (r^  are  the  Pauli-matrices  acting  on  the  Nambu-spinors): 

Hbcs  =  T,i^k  -  /^)*Ir3*fc  +  E  ^fc(9)*I+9/2n*fc-,/2  (10.1) 

k  k,q 

where  k  =  {p+p')/2  is  the  Fourier  transform  of  the  relative  coordinate  r  =  xi~X2 
of  two  interacting  electrons  and  q  =  p'  -  p  is  the  Fourier  transform  of  the  center 
of  mass  R  =  {xi  +  X2)/2.  Without  loss  of  generality  for  my  purposes  I  have 
assumed  a  real  order  parameter  Ak{q)  in  the  above  hamiltonian. 

Since  the  hamiltonian  is  bilinear,  the  corresponding  Green  function  can  be 
evaluated  and  has  the  simple  form  (for  Matsubara  frequencies  a;„  =  nT{2n  +  1) 
at  a  given  temperature  T;  I  drop  the  index  n  whenever  possible) 


(10.2) 
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Observe  that  due  to  the  Pauli  matrices  the  Matsubara  Green  function  is  itself  a 
2x2  matrix.  I  denote  this  by  the  hat.  The  theory  becomes  self-consistent  upon 
imposition  of  the  gap  equation,  i.e.  the  requirement,  that  the  order  parameter 
is  itself  an  expectation  value  of  a  pair  of  electrons.  In  general,  the  gap  equation 
reads 


The  tracing  of  the  product  ^Go  ensures  that  we  pick  the  correct  component  of 
the  matrix  Green  function.  Observe  that  the  right  hand  side  depends  implicitly 
on  Afc  via  the  Green  function. 

For  a  bulk  superconductor  without  impurities  the  order  parameter  is  suppos- 
edly uniform  in  space.  This  means  that  only  the  g  =  0  component  of  the  momen- 
tum dependent  order  parameter  is  nonzero,  Ak{q)  =  Afc(g  =  0)  =  A^.  As  a  re- 
sult, the  hamiltonian  10.1  can  be  diagonalized  (via  a  Bogoliubov-transformation) 
and  the  excitation  spectrum  of  quasiparticles  can  be  found.  It  is 


For  a  clean  isotropic  s-wave  order  parameter  A^  =  A^  the  spectrum  Ek  has  a  gap 
of  Ao  about  the  Fermi  level.  Therefore,  Ak  is  often  named  the  gap.  It  is  important 
to  distinguish  between  the  maximum  of  the  order  parameter  Afc  and  the  (maybe 
nonexistent)  gap  in  the  excitation  spectrum  for  all  anisotropic  superconductors. 
Although  related,  they  are  only  identical  for  the  isotropic  s-wave  SC's. 


(10.3) 


(10.4) 
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The  distinction  of  s-,  extended  s-,  p-  (e.g.  superfluid  ^He),  and  d-wave  order 
parameters  comes  solely  from  the  pairing  potential  Vk^k'  as  can  be  seen  from  the 
gap  equation  10.3.  For  simplicity,  I  use  a  pairing  interaction  which  factorizes  in 
the  momentum  variables  k,k': 

Vk,k'  =  T.^i^Kk)Mk')  (10.5) 

The  index  /,  /  =  0 . . .  oo  indicates  the  angular  momentum  "quantum  number"  of 
the  pair  wave  function,  whereas  $/(A;)  is  a  corresponding  normalized  basis  func- 
tion. The  pairing  interactions  Vi  are  taken  to  be  independent  of  the  magnitude 
of  k  for  energies  within  the  interval  [-uJo,uJo]  about  the  Fermi  level,  and  zero 
otherwise.  This  is  another  approximation  of  the  weak  coupling  BCS-theory.  The 
selection  of  an  order  parameter  with  a  specific  symmetry  is  achieved  by  setting 
all  Vi  to  zero  but  the  desired  one.  For  example,  42_y2-wave  superconductivity 
is  achieved  by  taking  V2  =  V  as  the  only  nonzero  coupling.  As  a  convenient 
choice  for  the  basis  function  for  a  cylindrical  Fermi  surface  in  two  dimensions 
one  can  choose  ^2{k)  =  y2cos(20fc),  where  0^  is  the  angle  of  the  momentum  k 
to  the  X-axis.  Because  of  the  momentum  factorization  of  the  pairing  interaction 
the  /c-dependence  of  the  right  hand  side  of  the  gap  equation  can  be  pulled  out 
of  the  sums.  It  follows  that  the  order  parameter  =  A^^  will  have  the  same 
angle  dependence  as  the  chosen  basis  function.  This  holds  even  in  the  case  of  a 
nonuniform  gap  in  the  case  of  an  impure  superconductor. 

With  a  separable  pairing  interaction  of  the  form  above,  I  can  partially  perform 
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the  momentum  sum  leading  to  (for  the  d-wave  case,  Ajt  =  Acos2(j)k) 


l  =  -2TVNo7rY^ 


— — cos  20fc 


1 


(10.6) 


The  sum  over  Matsubara  frequencies  is  constraint  by  [— cjojo;,,],  since  the  pairing 
potential  vanishes  otherwise.  No  is  the  density  of  states  at  the  Fermi  level  of  the 
normal  state. 

One  obtains  from  the  Eq.  10.6  by  setting  =  0.  The  result  is  identical 
to  the  one  quoted  above,  with  V  =  V2.  The  temperature  dependence  of  Ajt 
is  also  identical  to  the  standard  s-wave  case;  however,  the  value  at  T  =  0  is 
Afc(T  =  0)  =  2.UTc^2{k)  for  d-wave,  rather  than  Ak{T  =  0)  =  l.76Tc  for  the 
s-wave  case.  The  density  of  states  follows  from  the  order  parameter  and  the 
explicit  form  of  the  Green  function  Eq.  10.2  since 


where  Gi  1  is  the  upper  left  element  of  the  matrix  Green  function  G.  For  a  clean 
rfi2_j,2-wave  or  extended  s-wave  order  parameter  A^(a;)  is  linear  at  low  frequencies 
as  shown  in  figure  8.2.  These  types  of  superconductors  are  therefore  gapless  even 
in  the  clean  limit.  The  linear  frequency  dependence  of  N{u)  is  reflected  in  low 
temperature  thermodynamic  properties  like  the  London  penetration  depth  A(T). 


N{lo)  =  — Im y2Gii{k,iu  ^  to  +  ie) 


(10.7) 


k 


10.2    Self-Consistent  T-Matrix  Approximation 


The  self-consistent  T-Matrix  approximation  (SCTMA)  is  actually  not  a  defi- 


66 


nite  approximation  but  rather  a  term  given  to  the  type  of  approximations  in  which 
the  T-matrix  (defined  below)  of  a  single  impurity  problem  is  made  self-consistent 
in  order  to  treat  low  concentrations  of  impurities  for  which  one  expects  the  single 
impurity  scattering  to  be  dominant  over  complicated  multi-impurity  scattering 
processes.  Its  validity  for  three  dimensional  systems  is  well  established.  For  the 
layered  high-Tg  SC's,  for  which  one  might  expect  two-dimensional  physics  to  be 
relevant,  the  applicability  of  the  SCTMA  has  been  challenged  [67]  for  the  cases 
with  order  parameter  nodes.  I  will  argue  in  the  next  chapter  that  for  realistic 
disorder  in  a  two-dimensional  d-wave  SC  the  SCTMA  gives  in  fact  qualitatively 
correct  results.  It  is  therefore  safe  to  use  the  SCTMA  as  long  as  one  bears  in 
mind  that  certain  constants  might  need  experimental  input  rather  than  being  de- 
termined by  the  theory  itself  from  microscopic  parameters.  However,  qualitative 
physics  like  the  question  of  power  laws  can  be  addressed  by  the  SCTMA. 

Consider  nonmagnetic  impurities  in  a  bulk  superconductor  described  by  the 
Green  function  Eq.  10.2.  Static  impurities  (potential  scatterers)  are  described 
by  the  hamiltonian 


where  U{x-  Rimp)  is  the  scattering  potential  of  a  impurity  located  at  Rimp,  and 
U{p,p')  is  its  Fourier  transform.  Due  to  the  impurity  the  full  Green  function 
G  will  be  nonuniform  and  therefore  depend  on  two  momenta  p,p'.  Since  the 
impurity  is  static  (it  cannot  absorb  or  emit  energy)  G  will  still  be  a  function  of  a 


J:'^Hp)U{p,p')t,^{p') 


(10.8) 
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single  frequency  uj.  If  the  impurity  is  local,  that  is  U{x-Rimp)  =  UoS{x  —  Rimp),  a 
Dirac  (5-function  in  space,  the  full  Green  function  still  has  nontrivial  momentum 
dependence,  however,  the  irreducible  self  energy  will  be  a  function  of  frequency 
alone  if  certain  multi-impurity  scattering  processes  are  neglected  (irreducible 
means  that  no  diagrams  are  allowed  that  separate  into  two  disconnected  pieces 
upon  cutting  a  single  electron  line,  see  below).  In  general,  the  self  energy  E  is 
defined  via  Dyson's  equation 


G{uJ,p,p')  =  Go{io,p)5py  +  Go{uj,p)'t{p,uj)G{uj,p,p') 


(10.9) 


In  practice,  the  self  energy  is  computed  via  a  perturbation  expansion.  The 
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Figure  10.1:  Lowest  order  irreducible  self  energy  diagrams  for  scattering  off  local 
impurities.  The  cross  stands  for  the  impurity,  dashed  lines  indicate  the  interac- 
tions and  solid  lines  represent  the  electrons,  a),  b)  and  c)  correspond  to  scattering 
off  a  single  impurity,  d)  and  e)  are  diagrams  corresponding  to  two  impurity  scat- 
tering. 
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lowest  order  irreducible  diagrams  are  shown  in  Fig.  10.1.  The  cross  stands 
for  the  impurity,  dashed  lines  indicate  the  interactions  and  solid  lines  represent 
the  quasiparticles.  The  first  diagram  is  a  constant  and  can  be  absorbed  in  the 
chemical  potential.  The  second  diagram  corresponds  to  the  second  order  Born 
approximation.  Such  an  approximation  can  be  only  justified  for  small  interaction 
potential  U.  For  U  of  the  order  of  the  Fermi  energy  higher  order  diagrams,  e.g. 
Fig.  10.1  c),  become  important  and  have  to  be  accounted  for. 

This  can  be  done  elegantly  with  the  introduction  of  the  T-matrix.  In  the 
special  case  of  a  single  (5-function  scatterer  the  T-matrix  has  a  simple  form  and 
can  be  determined  as  follows:  The  T-matrix  is  defined  by  the  equation 


where  U {p,  p')  =  U {p,  p')t^.  For  a  (5-function  potential  U (p,  p')  =  [/„  is  a  constant. 
It  is  intuitively  clear  (and  can  be  shown  rigorously)  that  the  T-matrix  can  have 
no  momentum  dependence,  f{uj,p,p')  =  f{u).  The  momentum  sum  over  p"  can 
therefore  be  performed  and  leads  to 


,p,p')  =  U{p,p')  +  EU{p,p")Go{ico,p")f{u;,p",p') 


(10.10) 


f{u)  =  U0T3  +  UoTzgo{iio)f{u) 


(10.11) 


where  go{iuj)  =  J2pGo{iuJ,p)  is  the  momentum  integrated  Green  function  of  the 
bulk  SC.  For  a  general  SC,  the  solution  is 


f{u)  =  UoTsiTo  -  UoT^goiuj))  ^ 


(10.12) 
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with  To  the  2x2  identity  matrix.  For  a  d-wave  SC  the  angular  integral  of  the  off- 
diagonal  component  in  go{Lo)  vanishes.  The  T-matrix  then  has  the  particularly 


simple  form 


^^^^  -  ~T^mur '  ^''-''^ 


that  is,  there  is  no  off-diagonal  component  in  the  T-matrix  itself.  In  terms  of  the 
T-matrix,  the  full  Green  function  is  given  by 

G{ico,p,p')  =  Go{iuj,p)Sp^p'  +  Go{iuj,p)f{u)Go{iuj,p')  (10.14) 

As  advertised,  the  Green  function  still  strongly  depends  on  the  momenta  although 
neither  the  self  energy  nor  the  T-matrix  do  so. 

Rather  than  a  single  impurity  I  am  interested  in  a  low  concentration  of  im- 
purities in  a  bulk  sample.  The  concentration  must  be  low  enough  so  that  multi- 
impurity  scattering  processes  are  unimportant  compared  to  scattering  off  a  single 
impurity  site.  One  does  not  expect  the  actual  configuration  of  impurities  to  be 
of  any  importance  in  a  macroscopic  sample.  Therefore,  I  have  to  average  e.g. 
the  self  energy  diagrams  over  all  possible  configurations  of  impurities.  Details  of 
this  straightforward  procedure  can  be  found  e.g.  in  the  article  of  Ambegoakar  in 
Parks'  book  on  superconductivity  [76].  The  basic  result  is  that  the  disorder  av- 
eraged self  energy  and  the  T-matrix  are  simply  related  (for  (J-function  scatterers 
only)  by 

T,{lo)  =  niT{u)  +  crossed  diagrams  ,  (10.15) 
Ui  being  the  impurity  concentration.  In  a  three  dimensional  system,  due  to  phase 
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space  constraints  the  crossed  diagrams  obtain  a  factor  l/{kpl)  compared  to  the 
uncrossed  diagrams,  with  kp  the  Fermi  wave  number  and  /  the  elastic  mean  free 
path.  Since  typically  for  dilute  concentrations  of  impurities  kpl  »  1  the  crossed 
diagrams  can  be  neglected.  In  a  two  dimensional  system  there  is  no  such  a  priori 
suppression  of  the  crossed  diagrams.  However,  as  I  will  show  in  the  next  chapter, 
it  still  makes  sense  to  neglect  the  crossed  diagrams,  since  the  qualitatively  correct 
behavior  is  reproduced  by  this  approximation. 

The  crucial  point  in  Eq.  10.15  is  that  the  T-matrix  has  the  structural  form  of 
a  single  impurity  T-matrix,  like  the  one  of  Eq.  10.13.  However,  the  Qo  in  it  is  not 
the  one  obtained  from  the  bare  (noninteracting)  Green  function  but  rather  from 
the  Green  function  with  self  energy  E(a;).  This  means  that  Eq.  10.15  is  actually 
a  self-consistency  equation  for  the  disorder  averaged  (momentum  independent) 
self  energy  E(a;). 

The  full  solution  of  Eq.  10.15  is  usually  only  numerically  accessible,  although 
important  analytical  results  can  be  obtained  in  certain  parameter  regimes,  e.g. 
for  small  frequency  cj.  Having  obtained  the  self  energy,  one  can  determine  the 
DOS  and  from  there  the  thermodynamics  of  the  system.  Within  linear  response, 
one  can  also  obtain  information  about  the  electronic  contributions  to  transport 
properties  like  conductivity  and  thermal  conductivity. 


CHAPTER  11 
VALIDITY  OF  T-MATRIX  APPROXIMATION  IN  2D 

11.1    General  Problem 

Recently,  Nersesyan  et  al.  [67]  have  questioned  the  validity  of  the  SCTMA 
when  applied  to  a  strictly  two-dimensional  disordered  d-wave  SC.  They  pointed 
out  that  in  2D  no  small  factors  like  l/{kp'l)  distinguish  between  uncrossed  and 
crossed  diagrams.  In  fact,  straightforward  computation  shows  that  the  lowest 
order  crossed  diagram  Fig. 10.1  e)  (for  outer  momentum  A;  =  0)  is  of  the  same 
order  as  the  corresponding  uncrossed  diagram  Fig.  10.1  d)  that  is  accounted  for 
within  the  SCTMA,  namely  ~  to"^  log^(A/a;)  (A  being  the  maximum  of  the  order 
parameter).  This  is  only  true  for  outer  momentum  k  —  0,  since  finite  k  cut  off  the 
logarithmic  'divergence'  at  low  frequencies  for  the  crossed  diagram.  On  the  other 
hand  the  uncrossed  diagram  is  momentum  independent  since  the  integrations  over 
the  inner  momenta  qi,q2  are  completely  independent  of  k.  Although  this  caveat 
shows  that  even  in  2D  crossed  and  uncrossed  diagrams  are  not  on  equal  footing, 
it  is  clear  that  one  can  not  neglect  the  uncrossed  diagrams  without  establishing  at 
least  quaUtatively  in  an  independent  way  that  such  an  approximation  is  justified. 

Nersesyan  et  al.  avoided  perturbation  theory  by  applying  bosonization  tech- 
niques together  with  the  replica  trick.  They  found  a  power  law  DOS  N{u)  ~  \uj\°', 
a  ~  1/7,  for  sufficiently  small  frequency  u  and  weak  disorder,  rather  than  the 
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analytic  behavior  N{uj)  ~  const  +  aco'^  expected  within  SCTMA.  Their  calcula- 
tion supports  their  claim  that  the  uncrossed  diagrams  should  not  be  neglected, 
and  that  the  SCTMA  breaks  down  for  a  two-dimensional  d-wave  SC.  Although 
the  physical  systems  in  question  are  in  reality  highly  anisotropic  3D  systems,  the 
possibility  of  a  2D-3D  crossover  at  low  temperatures  could  conceivably  invali- 
date some  of  the  results  of  the  usual  "dirty  d-wave"  approach  using  the  SCTMA. 
This  would  render  the  description  of  the  low-temperature  transport  properties  of 
the  cuprate  superconductors  considerably  more  complicated  even  were  the  order 
parameter  of  the  simple  2D  dx2_y2  form  usually  assumed. 

In  contrast,  I  show  [56]  that  for  certain  types  of  disorder,  exact  results  can  be 
obtained  for  the  DOS  of  strictly  2D  disordered  superconductors.  I  show  that  for 
any  disorder  diagonal  in  position  and  particle-hole  space,  the  DOS  of  a  classic 
isotropic  s-wave  superconductor  has  a  rigorous  threshold  at  the  (unrenormalized) 
gap  edge  A,  as  expected  from  Anderson's  theorem  [53].  Within  the  same  general 
method,  I  show  that  the  residual  DOS  N{0)  of  a  superconductor  with  nodes  (e.g. 
d-  or  extended  s-wave)  is  nonzero  for  arbitrarily  small  disorder.  These  findings 
are  in  disagreement  with  Ref.  [67]  but  do  agree  qualitatively  with  the  results 
of  the  SCTMA  for  that  system  and  dimension.  Therefore,  I  conclude  that  the 
SCTMA  is  a  valid  method  to  qualitatively  describe  the  considered  system  and 
may  be  used  to  obtain  information  on  transport  properties. 

Nersesyan  et  al.  also  argued  that  a  non-zero  DOS  at  a;  =  0,  a  quantity  indi- 
cating spontaneous  symmetry  breaking,  may  not  occur  in  a  2D  system  because  of 
the  Mermin- Wagner  theorem  [77].  I  believe  that  the  DOS  in  a  disordered  system 
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is  not  an  order  parameter  which  belongs  to  the  class  of  order  parameters  covered 
by  the  Mermin-Wagner  theorem.  This  is  supported  by  the  fact  that  a  non-zero 
DOS  occurs  also  in  other  tight-binding  models  (e.g.,  model  for  two-dimensional 
Anderson  localization  [78]),  which  are  described  by  a  field  theory  with  continuous 
symmetry. 

11.2    Superconductor  on  a  Lattice  with  On-Site  Disorder 

The  method  of  calculating  exactly  the  DOS  of  a  superconductor  for  certain 
types  of  disorder  is  motivated  by  the  analysis  of  Dirac  fermions  in  2D  [79].  The 
BCS  Hamiltonian  in  first  quantized  form  is  given  by  (nV2m  =  1) 

H^{-V'-fi)T3  +  An  .  (11.1) 

It  describes  quasiparticles  in  the  presence  of  the  spin  singlet  order  parameter 
A.  As  before,  the  Tj  are  the  Pauli  matrices  in  particle-hole  (Nambu)  space.  The 
disorder  is  modeled  by  taking  the  chemical  potential  =  /Xj;  as  a  random  variable 
distributed  according  to  a  probability  distribution  P(/ii). 

I  consider  a  2D  square  lattice  spanned  by  the  unit  vectors  ii  and  62-  The 
kinetic  energy  operator  -V^  is  defined  by  its  action  on  a  wave  function  ^(x). 

VH{x)  =  '^{x  +  2ei)  +  ^(x  -  2ei)  +  <b{x  +  262)  +  *(x  -  262)  .  (11.2) 

This  definition  involves  displacements  of  two  lattice  sites  rather  than  one,  as 
would  be  the  case  in  the  simplest  tight-binding  representation  of  the  lattice  kinetic 
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energy.  The  reason  for  this  choice  is  technical  convenience  at  a  later  part  of  the 
calculation.  For  a  system  of  fermions  in  the  thermodynamic  limit,  the  kinetic 
energy  will  have  a  band  representation  quite  similar  to  the  usual  tight-binding 
form.  In  particular,  there  will  be  no  distinguishing  features  of  the  band  structure 
near  the  Fermi  level.  The  definition  obeys,  of  course,  the  same  global  continuous 
symmetries  discussed  for  the  model  in  Ref.  [67]. 

The  bilocal  lattice  operator  A  =  A^:  is  taken  to  act  as  a  c-number  in  the 
isotropic  s-wave  case, 

A^{x)  =  A^(a;),  (11.3) 
whereas  to  study  extended  (nonlocal)  pairing  I  define 

A'd^bix)  =  A^ [^(x  +  ei)  +  ^{x  -  ei)  ±  ^{x  +  63)  ±  "^{x  -  62)] .  (11.4) 

These  are  the  standard  representations  of  the  corresponding  order  parameters  on 
a  square  lattice. 

I  consider  the  single-particle  Matsubara  Green  function  defined  as  G{iu)  = 
{iuTo  -  H)~^  (suppressing  the  index  on  the  Matsubara  frequencies  a;„).  I  am 
interested  in  calculating  the  DOS 

A'' (a;)  =  -^lm{Gi-i{x,x,iuj      u  +  ie))  (H-S) 

where  (...)  denotes  the  disorder  average.  The  major  problem  is  how  to  perform 
this  disorder  average  over  the  probability  measure  P{n^)dn^  of  the  random  vari- 
able /ij.  Exact  results  for  the  disorder-averaged  Green  function  in  noninteracting 
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systems  can  frequently  be  obtained  for  Lorentzian  disorder,  by  exploiting  the 
simple  pole  structure  of  in  the  complex  fi^  plane. 

P{li,)di^^  =  '  \'      2  ^Mx ,  (11.6) 

Ho  is  the  chemical  potential  of  the  averaged  system.  For  convenience,  I  set  /io  =  0. 
The  averaged  Green  function 

(GiiLo))  =  f  X{d^Ji,P{^i,)G{^uJ■^i,)  (11.7) 

•'  X 

can  be  easily  evaluated  if  G  can  be  shown  to  be  analytic  in  either  the  upper  or 
lower  /i-half  plane. 

In  a  superconductor,  the  Green  function  depends  on  the  random  variable 
via  ±  i'jJ-  This  is  a  consequence  of  the  particle-hole  structure,  i.e.  the 
different  Pauli  matrices  multiplying  cj  and  Therefore,  the  averaging  of  G  with 
respect  to  Lorentzian  disorder  is  not  trivially  possible.  However,  it  is  possible  to 
reformulate  the  problem  so  that  G  is  a  sum  of  terms  each  of  which  are  analytic 
in  either  the  upper  or  the  lower  complex  ix^  plane.  This  allows  then  to  perform 
the  averaging  of  the  Green  function  for  Lorentzian  disorder. 

11.3    Isotropic  S-Wave  Superconductor 

Consider  first  a  homogeneous  and  isotropic  s-wave  order  parameter.  The 
Matsubara  Green  function  may  be  written  G{iijj)  -  -{iutq  +  H){uj'^  +  j 
note  that  H'^  =  (- V^-;/)Vo+ AVq  since  in  the  isotropic  s-wave  case,  (-V2-/i)r3 
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anticommutes  with  Ati  even  for  random  /i^.  This  is  due  to  the  locality  of  the 
order  parameter  A  in  this  case. 

The  expression  H'^+lu'^tq  is  proportional  to  the  unit  matrix;  as  a  consequence, 
the  Green  function  can  be  written  in  the  simple  form 


G{iuj)  —  — 


2iVA 


2  ^^2 


1 


-W^  -  H  -  i^/^^nJ^     - V2  - // +  zx/A2  +  a;2 


To 
(11.8) 


The  imaginary  part  of  this  expression  (after  analytic  continuation)  for  any  given 
configuration  of  impurities  is  vanishing  for  |a;|  <  A.  Therefore,  the  DOS  shows 
a  gap  of  size  A  independent  of  the  distribution  function  P(//).  Thus,  the  model 
reproduces  the  famous  Anderson  theorem  [53]  which  states  that  the  thermo- 
dynamics of  an  isotropic  s-wave  superconductor  are  not  affected  by  diagonal, 
nonmagnetic  disorder.  The  situation  is  different  if  the  order  parameter  itself  is 
random  [80,81].  In  that  case  all  quasiparticle  states  are  broadened  and  the  DOS 
is  finite  even  for  the  isotropic  s-wave  SC. 

11.4    D-  and  Extended-S  Symmetry  Superconductors. 


My  main  concern  is  with  the  d-wave  and  extended-s  "  bond"  order  parameters 
Ad  defined  above.  The  corresponding  pure  systems  in  momentum  space  fulfill 
the  condition  Y.k  =  0,  so  that  nonmagnetic  disorder  must  cause  significant 
pair  breaking  [54].  The  behavior  of  the  imaginary  part  of  the  Green  function  can 
be  studied  using  a  method  analogous  to  that  used  for  the  s-wave  case.  However, 
because  of  its  nonlocal  nature  the  order  parameter  Adri  does  not  anticommute 
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with  (— —  fi)T3  anymore  if  is  random.  As  a  consequence,  is  no  longer 
proportional  to  the  unit  matrix  tq,  which  forbids  a  simple  separation  of  poles 
of  the  Green  function  as  in  the  s-wave  case,  Eq.  11.8.  A  different  type  of 
transformation  is  required.  I  introduce  a  diagonal  matrix  (or  staggered  field) 
Dx,x'  =  {-^y^'^'^'^^x,x'  (note  D'^  is  the  unit  matrix).  Now  I  may  write 

=  HDtIDH  =  [(-V^  -  /x)Dro  -  i^'-^ Dt2][D{-V''  -  h)to  +  zDA^ra]  (11.9) 

D  commutes  with  -V^  and  /i  as  defined  above,  because  -V^  involves  only  next 
nearest  neighbor  sites  and  /i  is  local.  However,  D  anticommutes  with  the  order 
parameter  ,  since  the  nonlocal  order  parameter  involves  nearest  neighbor  sites. 
This  yields  simply       =  H"^,  with 

H  =  (-V^  -  ij,)Dto  -  iAWt2  .  (11.10) 

Therefore,  the  quantity  +  cjVq  =  (H  +  iuto){H  -  iuTo)  can  be  used  to 
write 

i{iLOTo  +  H)  f  1  1  \ 

GM  =  ^^  -[Yr  ■■  7v — ■   11-11 

Note  that  both  E  and  E  appear  in  this  expression,  but  E  only  in  the  numerator. 

Defining  =  fJ.xDx,x,  I  now  note  that  for  a;  >  0  and  lm{zx)  >  0  the  matrix 
icjTo  -  E  is  non-singular  (i.e.,  det{iujTo  ~  E)  ^  0).  Therefore,  the  transformed 
Green's  function  {iutq  -  Ey^  can  be  expanded  as  a  Taylor  series  with  nonzero 
radius  of  convergence  around  any  Zx  in  the  upper  half  z-  plane,  and  is  conse- 
quently analytic  there.  Correspondingly,  {iutq  +  E)'^  is  analytic  in  the  lower 
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half  z—  plane.  Due  to  these  analytic  properties  and  using  P{z^  =  I  can 

now  straightforwardly  perform  the  disorder  integration.  Care  has  to  be  taken 
when  evaluating  the  numerator  involving  i/,  since  it  involves  rather  than  z^.. 
It  turns  out,  though,  that  all  terms  involving  the  matrix  D  vanish. 

The  resulting  disorder  averaged  Matsubara  Green  function  is  translational 
invariant.  Performing  a  spatial  Fourier  transform  I  replace  —  by  ^  =  -  /xq 
and  obtain 

This  is  the  Matsubara  Green  function  of  the  pure  system  with  the  frequency 
iu)  shifted  by  the  disorder  parameter,  iu)  ico  +  ij.  In  contrast,  for  the  local 
(isotropic)  s-wave  order  parameter  discussed  before,  the  average  over  a  Lorentzian 
distribution  in  Eq.  11.8  implies  a  shift  iy/A'^  +  a;^      iy/A^  +  uP'  +  Z7. 

To  obtain  the  DOS  for  the  d-wave  case  I  approximate  the  sum  over  the 
momenta  k  in  standard  fashion  as  No  Jq^  ^  where  No  is  the  density  of 

states  of  the  normal  metal  at  the  Fermi  level.  I  also  approximated  the  tetragonal 
Fermi  surface  of  a  square  lattice  by  a  circle.  The  result  is 

where  the  d-wave  order  parameter  is  approximated  by  Ad((?!»)  =  cos (20).  At 
u  =  0,  this  leads  to  A^(0)  =  A^o^ln(4Ad/7)  for  7  <<  A^.  Thus,  the  density 
of  states  is  nonzero  at  the  Fermi  level  for  arbitrarily  small  values  of  the  disorder 
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parameter  7.  If  I  expand  the  integral  for  small  values  of  I  find  that  N{u;)  rises 
as  w^. 

For  more  general  continuous  distributions  P{fj,)dfi  the  averaged  DOS  can 
be  estimated  using  again  the  analytic  structure  of  G.  Applying  the  ideas  of 
Ref.  [80],  one  can  derive  a  lower  bound  by  a  decomposition  of  the  lattice  into 
finite  sub-blocks.  The  average  DOS  on  an  isolated  sub-block  can  be  estimated 
easily.  Moreover,  the  contribution  of  the  connection  between  the  sub-blocks  to 
the  average  DOS  can  also  be  estimated.  A  combination  of  both  contributions 
leads  to  (iV(0))  >  cj  min_^j<^<^j  P(/x),  where  ci  and  //i  are  distribution  depen- 
dent positive  constants.  In  particular,  fii  must  be  chosen  such  that  the  spectrum 
of  H{no  =  0)  =  -V^r3  +  AdTi  is  inside  the  interval  [-1^1,^1].  For  all  unbounded 
distributions,  like  the  Gaussian  distribution  used  in  Ref.  [67],  as  well  as  compact 
distributions  with  sufficiently  large  support  this  estimate  leads  to  a  nonzero  DOS 
at  the  Fermi  level. 

11-5    Consequences  and  Comparison  to  Other  Methods. 

The  major  result  in  the  d-wave  (extended  s-wave)  case  with  Lorentzian  dis- 
order is  the  presence  of  a  finite  purely  imaginary  self-energy  Eq  =  -^7ro  due  to 
nonmagnetic  disorder.  This  leads  to  a  nonzero  DOS  at  the  Fermi  level,  in  qual- 
itative agreement  with  standard  theories  based  on  the  SCTMA  [57,  58]  as  well 
as  with  exact  diagonalization  studies  in  2D  [82].  In  contrast  to  such  theories  the 
above  self  energy  has  no  dependence  on  A2,  i.e.  it  is  the  same  as  in  the  normal 
state.  In  Fig.  11.11  show  a  comparison  of  the  self  energies  of  my  theory  and  the 
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weak  (Born)-  and  unitary-scattering  limits  of  the  SCTMA. 
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Figure  11.1:  Imaginary  part  of  the  self  energy  vs.  frequency.  For  Lorentzian 
disorder  (solid  line)  the  self  energy  is  constant  27.  The  self  energy  of  the  self- 
consistent  T-matrix  approximation  in  the  unitary  scattering  limit  (dashed-dotted 
line)  behaves  oc  (A^)^/^  at  zero  frequency.  For  Born  scattering  (dashed  line)  the 
value  at  zero  frequency  is  nonzero,  but  exponentially  small.  I  have  adjusted  the 
impurity  concentration  to  obtain  equal  normal  state  self  energies  for  the  T-matrix 
results. 

A  drawback  of  the  model  with  Lorentzian  disorder  is  that  impurity  concentra- 
tion does  not  appear  explicitly  in  the  theory.  Whereas  in  the  T-matrix  approach 
I  have  with  the  impurity  concentration  and  the  scattering  strength  (or  phase 
shift)  two  parameters  associated  with  disorder,  in  the  present  model  I  have  only 
7,  the  width  of  the  Lorentzian.  A  way  of  making  a  connection  is  by  comparing 
the  variance  of  the  Lorentzian  distribution  (7)  and  the  variance  of  the  distribu- 
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tion  underlying  the  T-matrix  approximation,  which  is  a  bimodal  distribution  of 
a  chemical  potential  ^  =  hq  with  probability  \  —  6  {6  being  the  dimensionless  im- 
purity concentration)  and  ix  —  ijlq-\-V  with  probability  5  {V  being  the  scattering 
potential).  The  variance  Var^  of  this  distribution  is  determined  by 

V^rl  =  {^?)-{^,f  =  V\5-5^).  (11.14) 

For  small  concentrations  of  impurities,  6  «  1,  I  find  Var^  =  VS^^'^.  The  5^^'^ 
behavior  is  also  found  for  ImEo(a;  =  0)  in  the  T-matrix  approach  for  strong 
scattering.  Since  in  my  model  the  variance  of  the  distribution  is  also  the  imagi- 
nary part  of  the  self  energy,  this  suggests  that  my  model  is  closer  to  the  strong 
scattering  limit  of  the  SCTMA  than  the  Born  limit. 

Finally,  I  comment  on  the  discrepancies  between  my  result  and  the  calculation 
of  Nersesyan  et  al.,  who  found  a  power  law  for  the  averaged  DOS  with  Gaussian 
disorder. 

One  might  question  the  analysis  of  Nersesyan  et  al.  because  of  the  use  of  the 
replica  trick,  which  is  a  dangerous  procedure  in  a  number  of  models.  [83]  How- 
ever, Mudry  et  al.  [84]  have  obtained  identical  results  for  the  continuum  problem 
of  Dirac  fermions  in  the  presence  of  a  random  gauge  field  using  supersymmetry 
methods.  I  therefore  believe  that  the  crucial  difference  between  my  results  and 
those  of  Ref.  [67]  occurs  in  the  passage  to  the  continuum  and  concomitant  map- 
ping of  the  site  disorder  in  the  original  problem  onto  the  random  gauge  field. 
Only  in  the  continuum  case  is  there  a  direct  analogy  between  disorder  in  the 
chemical  potential  and  a  gauge  field;  on  the  lattice,  gauge  fields  and  chemical 


82 


potential  terms  enter  quite  differently.  First,  chemical  potential  terms  are  local 
while  gauge  fields  are  defined  on  bonds.  Furthermore,  chemical  potential  disor- 
der enters  linearly  in  the  Hamiltonian  while  gauge  fields  enter  through  the  Peierls 
prescription  as  a  phase  in  the  exponential  multiplying  the  kinetic  energy. 

Disorder  of  the  gauge  field  type  is  furthermore  nongeneric  even  in  the  con- 
tinuum, as  discussed  by  Mudry  et  al.,  who  showed  that  the  critical  points  of 
the  system  with  random  gauge  field  are  unstable  with  respect  to  small  perturba- 
tions by  other  types  of  disorder.  I  expect  that  a  proper  mapping  of  the  lattice 
Dirac  fermion  or  d-wave  superconductor  problems  to  continuum  models  will  in- 
evitably generate  disorder  other  than  random  gauge  fields.  Therefore,  I  believe 
that  my  result  of  a  finite  DOS  at  the  Fermi  level  is  the  generic  case  for  a  d-wave 
superconductor  in  two  dimensions. 

In  summary,  this  calculation  suggests  that  the  standard  T-matrix  approach 
to  disordered  d-wave  superconductors  is  qualitatively  sufficient.  I  doubt  that  the 
result  by  Nersesyan  et  al.,  who  found  a  power  law  for  the  averaged  DOS  with 
Gaussian  disorder,  is  of  any  relevance  for  the  system  under  consideration. 


CHAPTER  12 
LOCAL  ORDER  PARAMETER  PERTURBATIONS 


Order  parameter  perturbations  in  the  vicinity  of  nonmagnetic  impurities  are 
present  even  in  the  classic  s-wave  SC's.  Fetter  [85]  has  calculated  the  pertur- 
bations about  a  single  impurity  within  a  continuum  model  in  three  dimensions, 
and  found  that  they  have  the  following  general  behavior:  a)  They  are  oscillatory, 
in  analogy  to  Friedel  oscillations  of  interacting  metals,  and  b)  they  decay  like 
a  power  law  within  a  distance  of  the  coherence  length  ^  and  exponentially  for 
distances  larger  than  ^.  A  sketch  of  this  behavior  is  shown  in  Fig.  12.1.  Later, 

1^ 


Figure  12.1:  Sketch  of  the  order  parameter  perturbations  around  a  nonmagnetic 
impurity.  The  perturbations  are  oscillatory  on  atomic  lengths  scales  a  and  decay 
exponentially  beyond  the  coherence  length  ^.  is  the  maximum  of  the  bulk 
order  parameter. 

Shiba  [86],  Rusinov  [87]  and  Schlottmann  [88]  have  worked  on  various  aspects  of 
the  problem  of  magnetic  and  nonmagnetic  impurities  in  an  s-wave  SC.  Recently, 
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numerical  work  has  been  performed  on  a  model  of  nonmagnetic  impurities  on  a 
two-dimensional  lattice  for  both  s-wave  and  d-wave  SC's  [82,89]. 

Despite  the  above  general  similarities,  there  are  important  differences  between 
the  order  parameter  perturbations  of  classic  s-wave  SC's  and  the  high- Tc  d-wave 
SC's.  These  are: 

•  The  order  parameter  perturbations  are  anisotropic,  reflecting  the  d-wave 
nature  of  the  underlying  SC. 

•  The  coherence  length  ^  is  much  shorter  than  in  high  Tc  SC's  than  in  classical 
SC's,  ^  ~  10  —  I5A,  which  is  a  few  atomic  distances. 

•  The  small  momentum  transfer  component  of  the  order  parameter  perturba- 
tions can  be  large  in  d-wave  SC's,  in  stark  contrast  to  the  isotropic  s-wave 
SC's  where  it  is  vanishing. 

The  first  point  is  not  surprising  and  of  minor  importance.  The  short  range  of 
the  perturbations  is  helpful  since  it  allows  us  to  treat  them  as  local  perturbations 
around  the  impurities  with  no  direct  interference  with  perturbations  due  to  other 
impurities  (in  the  dilute  limit).  It  is  the  third  point,  however,  which  makes  all  the 
diflference.  The  fact  that  the  d-wave  order  parameter  has  zero  average  when 
integrated  over  the  Fermi  surface  leads  to  a  qualitative  different  behavior  of  the 
d-wave  SC  even  in  the  case  of  a  single  impurity.  For  an  s-wave  SC  the  spatial 
integral  of  the  order  parameter  perturbations,  i.e.  the  q  =  0-component,  average 
out  to  zero,  5l^{q  =  0)  =  0,  in  agreement  with  Anderson's  theorem.  However, 
for  a  d-wave  SC  the  order  parameter  perturbations  have  a  nonzero  average. 
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Furthermore,  and  somewhat  counterintuitively,  6A{q  =  0)  is  not  bounded  by 
the  bulk  gap  A  [90]  as  I  will  show  below.  It  is  therefore  important  to  determine 
these  perturbations  self-consistently  via  the  gap  equation.  In  principle,  the  gap 
equation  asks  for  such  self-consistency  in  both  s-wave  and  d-wave  cases,  but  the 
vanishing  of  the  6A{q  =  0)  without  self-consistency  implies  that  the  corrections 
are  negligible  if  self-consistency  is  taken  into  account  for  the  s-wave  case.  I 
discuss  this  in  more  detail  below. 

12.1    Single  Impurity  Scattering,  Self-Consistent  to  First  Order 
For  a  bulk  d-wave  SC  the  (matrix)-Green  function  reads  (Eq.  10.2): 

Go{k,lLj)  =  {iLU-^T3-A^Ti)-\  (12.1) 

where  I  assume  an  order  parameter  of  the  form  A<^  =  A  cos  2(/>.  A  helpful  quantity 
is  the  momentum  integrated  Green  function  Qo  given  by 

So(-)  =  ET.-I^Go(t.^))  =  /'I^,.,  (12.2) 

k        2  Jo     27r(A;- (2Cj)2)l/2 

Note  that  there  is  no  off-diagonal  part  (a  ti  )  due  to  the  d-wave  symmetry.  For 
scattering  off  a  single  (^-function  impurity  with  strength  I  derived  in  chapter 
10  the  T-matrix,  Eq.  10.13 
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In  general,  the  (weak  coupling)  gap  equation  reads  (see  chapter  10,  Eq.  10.3) 

Ak{q)  =  -r  E  E  ^fc,fc'Tr  {riG{k'  +  q/2,  k'  -  q/2,  iu;)}  .  (12.4) 

1^  k' 

with  G  now  given  hy  G  =  Go  +  Go  +  GoTGo  as  in  Eq.  10.14.  Assuming  the 
pairing  potential  is  separable  and  has  d-wave  symmetry  [Vk^k'  =  V<^2{k)^2{k')) 
I  can  ask  for  the  change  of  the  gap  SAk{q)  due  to  the  impurity  scattering. 

5Ak{q)  =  -V^,{k)T  E  E  Mk'm  { "^Goik'  +  q/2,  iio)f{uj)Go{k'  -  q/2,  lu)  ]  . 

UJ      k'  ^   ^  ' 

(12.5) 

This  equation  determines  the  order  parameter  perturbations  to  zeroth  order,  in 
the  sense  that,  although  I  use  the  gap  equation,  I  assumed  an  unperturbed  order 
parameter  on  the  right  side  of  Eq.  12.5.  Clearly,  this  is  inconsistent. 

Shiba  and  Rusinov  [86, 87]  first  treated  the  perturbations  self-consistently  to 
first  order  (for  an  s-wave  SC,  for  d-wave  see  [90]  and  below).  They  noted  that, 
to  first  order,  it  is  enough  to  add  a  term  6Tk'{q,uj) 

STk>{q,  uj)  =  Go{k'  +  q/2,  iuj)5Ak'{q)Go{k'  -  q/2,  lu)  (12.6) 

to  the  r.h.s.  of  Eq.  10.13.  This  is  basically  the  first  term  in  a  perturbation  series 
for  the  T-matrix  with  a  scattering  potential  dAk'{q).  Now  both  the  left  and  right 
side  of  the  gap  equation  are  linear  in  SAk{q).  Due  to  linearity  one  can  separate 
the  terms  oc  SAk{q)  from  the  terms  not  containing  dAk{q).  For  brevity,  I  consider 
here  only  the  case  g  =  0,  leaving  the  g-dependent  expression  to  the  Appendix  C. 
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Since  5Ak{q  =  0)  =  cos  20^  the  A;-  dependence  drops  out  of  Eq.  12.5  with  the 
modified  T-matrix.  The  sum  over  momenta  k'  is  rewritten  in  standard  fashion, 
Efe  f  =  ^o!  <^?/  f^)  with  No  the  normal  density  of  states  at  the  Fermi 
level.  Only  few  terms  survive  the  tracing  procedure  and  the  integration  over  the 
energy  ^  =  e^/  —  //.  Explicitly,  the  equation  for  5 A  reads 


6 A  {\/V  +  2TY.Noj  ^cos'20  j 


(a;2  +  ^2  + 

The  ^-  integrals  are  easily  computed: 


2rE../gco.2,/^^-«^.  (...) 


r°°         1  TT  r°°        3  IT 


This  yields  for  Eq.  12.7 


a;2 


Adding  and  subtracting  in  the  numerator  of  the  integral  on  the  left  hand  side, 
I  obtain  a  term 

which  can  be  seen  by  comparing  with  the  gap  equation  10.6.   The  remaining 
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expressions  do  not  explicitly  depend  on  the  coupling  V.  The  final  result  is 


To  get  an  idea  about  the  order  of  magnitude  of  SA  I  consider  the  Ginzburg- 
Landau  (GL)  regime,  that  is,  temperatures  close  to  Tc,  so  that  the  order  param- 
eter is  small.  I  therefore  can  neglect  in  the  denominators  of  Eq.  12.11. 
This  renders  the  angular  integrals  trivial.  With  the  same  assumptions  I  find 
go  =  -iNoTT.  I  am  now  in  a  position  to  perform  the  Matsubara  sums.  The 
relevant  sums  are 

Tj:i/\^\'  =  c,/T\   T^l/a;^  =  c2/r,  (12.12) 

with  dimensionless  constants  C2,  C3  that  can  be  expressed  in  terms  of  Riemann 
C-functions.  Since  I  am  in  the  GL-regime  I  replace  the  temperature  T  by  Tc  and 
find 

[Nq-kUq)-^  C2 

^^=-'T;ty^-  (12.13) 

A  similar  result  was  obtained  by  Choi  and  Muzikar  [91]  in  a  different  context. 

12.2    Discussion  of  the  First  Order  Result 

This  result  has  been  deceptively  written  to  meet  the  expectation,  namely 
that  5A  is  some  negative  coefficient  times  the  bulk  gap  A.  However,  the  'coef- 
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ficient'  is  itself  inversely  proportional  to  A^.  Therefore,  SA  is  actually  inversely 
proportional  to  A.  Since  A  approaches  zero  as  T  ->  like  (1  -  T/TcY^"^,  5 A 
diverges  like  (1  -  T/Tc)~^^'^.  This  at  first  hand  counterintuitive  result  could  be 
understood  by  considering  that  the  g  =  0  component  of  6Ak{q)  is  nothing  but 
the  total  volume  integral  of  the  order  parameter  perturbations.  The  divergence 
oi6Ak{q  =  0)  therefore  signifies  that  the  order  parameter  perturbations  are  long 
ranged.  This  would  be  in  agreement  with  the  early  mentioned  general  behavior 
of  the  order  parameter  perturbations,  namely  that  their  range  (after  which  they 
decay  exponentially)  is  governed  by  the  coherence  length.  Since  the  coherence 
length  diverges  as  one  approaches  the  critical  point  so  does  the  range  of  the  order 
parameter  perturbations,  and  their  volume  integral  therefore  diverges. 

This  argument  seems  reasonable,  however,  it  neglects  that  the  approximation 
I  used  assumed  order  parameter  perturbations  to  be  small,  so  that  I  could  use 
the  standard  T-matrix  with  a  linear  correction  term.  A  diverging  SAk{q  =  0) 
invalidates  the  basic  assumption.  In  fact,  there  is  no  physical  reason  why  6Ak{q  = 
0)  should  diverge.  The  argument  above  shows  that  it  can  be  large.  Therefore, 
higher  order  calculations  are  indispensable,  at  least  in  the  unitary  scattering 
limit,  NonUo  »  1,  and  close  to  T^.  My  method  of  incorporating  all  orders  of 
6Ak{q  =  0)  (see  next  chapter)  leads  to  a  vanishing  SAk{q  =  0)  close  to  Tc,  rather 
than  diverging  one.  In  the  other  hand,  my  approximation  assumes  short  ranged 
order  parameter  perturbations,  an  assumption  that  is  clearly  violated  close  to  T^. 
At  the  moment,  it  is  unclear  to  me  whether  long  range  effects  or  higher  order 
corrections  dominate  close  to  Tc. 


90 


However,  I  am  not  really  interested  in  temperatures  close  to  Tc  but  the  effect 
of  the  order  parameter  perturbations  on  low  temperature  transport.  At  low 
temperatures,  the  linear  approximation  is  not  diverging,  but  still  SAk{q  =  0) 
turns  out  to  be  too  large  to  be  used  with  good  conscience.  I  therefore  treat 
SAk{q)  properly  to  all  orders.  This  can  be  achieved  by  considering  it  a  new, 
induced  scattering  term  in  the  hamiltonian  in  addition  to  the  standard,  direct 
scattering  term  of  the  nonmagnetic  impurities.  This  asks  for  a  completely  new 
determination  of  the  T-matrix.  In  the  next  chapter,  I  will  show  a  way  to  obtain 
an  analytical  result  for  the  T-matrix  that  incorporates  the  major  physics  of  the 
new  scattering  terms. 

To  conclude  this  chapter  let  me  stress  again  the  contrast  between  isotropic 
s-wave  and  d-wave  SC's  when  it  comes  to  SAk{q  =  0).  To  linear  order  one  finds 
that,  rather  than  diverging  in  the  GL-limit,  the  order  parameter  perturbations 
SAk{q  =  0)  =  0  for  all  temperatures.  This  is  because  for  an  s-wave  SC  there 
is  an  off-diagonal  term  (oc  tj)  in  the  momentum  integrated  Green  function  g^. 
This  term  turns  out  to  cancel  exactly  the  term  on  the  right  hand  side  of  the 
equivalent  of  Eq.  12.11.  Since  the  total  right  hand  side  of  the  equation  is  zero, 
so  must  be  SAk{q  =  0)  for  all  temperatures  T.  Within  my  method  of  calculating 
SAk{q  =  0)  to  all  orders  I  also  find  that  SAkiq  =  0)  =  0  is  a  solution  to  the 
nonlinear  equation,  though  other  solutions  might  be  possible.  I  therefore  believe 
that  for  an  s-wave  SC  SAk{q  =  0)  is  either  zero  or  too  small  to  play  any  relevant 
role. 


CHAPTER  13 


A  NEW  T-MATRIX  AND  CONSTITUTIVE  EQUATION  FOR  THE  ORDER 

PARAMETER  PERTURBATIONS 


13.1    General  Remarks 

It  was  shown  in  the  last  chapter  that  SAk{q)  is  not  a  priori  small  when  com- 
pared to  the  Fermi  energy  and  standard  scattering  terms.  It  is  therefore  necessary 
to  take  more  than  just  the  first  order  of  SAk{q)  into  account.  This  poses  a  serious 
problem,  since  SAk{q)  is  both  anisotropic  and  g-dependent  (nonlocal).  Especially 
the  ^-dependence  makes  the  problem  analytically  intractable,  since  dAk{q)  is  de- 
termined by  a  nonlinear  self-consistency  equation.  Future  work  might  use  either 
heavy  numerical  or  variational  techniques  to  improve  upon  this  calculation.  As 
a  first  step,  however,  I  try  to  obtain  approximate  results  by  exploiting  the  short 
range  of  the  fluctuations. 

The  anisotropy  (A;-dependence)  can  be  handled  as  soon  as  the  assumption 
of  locality  (no  ^-dependence)  is  made.  Assuming  cylindrical  symmetry  of  the 
Fermi  surface  I  can  perform  a  partial  wave  expansion  of  the  functions  on  the 
Fermi  surface,  that  is,  I  express  them  in  terms  of  the  coefficients  of  cos  Icj)  and 
sin/0.  I  assume  that  only  s-  and  d-wave  scattering  takes  place.  This  is  essentially 
an  assumption  about  the  pairing  potential  Fpy :  If  it  contains  only  s-  and  d-wave 
components  no  otlier  terms  can  be  induced.  In  principle,  the  technique  can  be 
extended  to  higher  partial  waves  (at  least  at  q  =  0). 
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The  procedure  described  first  computes  an  analytical  form  of  the  (now  mo- 
mentum dependent)  T-matrix  as  a  function  of  parameters  characterizing  diagonal 
(oc  Ta)  and  off-diagonal  (oc  6-  and  d-wave  scattering.  These  parameters  can 
then  be  determined  by  use  of  the  gap  equation.  This  can  only  be  done  numerically 
and  the  results  will  be  shown  in  the  next  chapter. 

13.2    Determination  of  the  T-matrix 

For  clarity  I  repeat  the  definitions  of  the  Green  function  and  its  momentum 
integral. 


Go{k,Lu)  =  {ioo  -  ^Ts  -  A^Ti) 


-1 


(13.1) 


(13.2) 


(13.3) 


For  a  bulk  d 


wave  SC  I  take       =  A  cos  2(j)  from  which  follows  that 


'2'^  d(f)  iu! 


(13.4) 


'o    27r  (A2  -  (za;)2)i/2 
dcf)  cos  20A  cos  2(f) 


01 


Jo 


0    27r(A2  -  (ia;)2)i/2 


(13.5) 


0  for  all  I  odd 


(13.6) 


I  also  have  =  The  higher  even  ^^'s  can  be  related  to  the  given  ones  and 
are  in  general  suppressed  by  a  factor  of  1//. 
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I  attempt  to  find  the  T-matrix  defined  by  the  equation  [92] 

Tpy{Lj)  =  U,,,,+       ^f/p,p,G'„((/>,a;)rp,y(u;)  (13.7) 

where  0  is  the  angle  of  pi  to  the  x-axis.  The  momenta  (p,p')  are  pinned  at  the 
Fermi  surface.  Up^p'  is  a  general  scattering  term  which  can  involve  diagonal  and 
off-diagonal  s-  and  d-wave  scattering.  I  expand  U,  Go  and  f  in  the  following 
way  (partial  wave  expansion): 

Up,p'  =  ipli'  cos  I4>p  cos  /Vp'  +  Ui  i,  sin  /0p  sin  /  Vp')  (13.8) 
G'o(0,  ^)  =  Y^  (jrn  COS  mt/)  (13.9) 

m 

%,p'  =  J2  {Tn,n'  COS  n(/)p  cos  n'(j)p'  +  f     sin ncpp  sin  nVp-)  (13.10) 

n,n' 

Here,  the  superscripts  e/o  denote  even  and  odd  terms  upon  reflection  on  the 
X-axis,  e.g.  0p      -cfjp.  Specifically,  I  choose  the  Ui^i'  =  lJi6i^i>  with 

Ul  =  Uo  =  UoTs  +  6,n  U°  =  0  (13.11) 

f>r  =  U",  =  ~5dn  (13.12) 

ijf  =  U^  =  Q  for  all  />1  ,  (13.13) 

leading  to 


Up^p'  -  UoT'i  +  6sTi  +  5d  (cos  (t)p  cos      -  sin  (j)p  sin  </)p')  ri . 


(13.14) 
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The  choice  for  Ui  is  motivated  by  considering  a  d-wave  scattering  term  like 
Sd  cos  2(j)k.  Since  k  =  {p  +  p')/2  I  have 

5d  cos  2(l)k  =  6d  (cos  (j)p  cos  (j)p'  -  sin  (/»p  sin  ^p/)  (13.15) 

for  forward  scattering  (g  =  |p  -  p'|  ->  0).  An  explicit  ^-dependence  is  neglected, 
since  at  the  moment  I  am  only  interested  at  — >  0. 

The  integral  (/)  in  Eq.  13.7  can  now  be  performed  at  the  cost  of  the  sum- 
mations over  the  partial  wave  indices  (which  run  over  all  integers).  The  result  is 
(suppressing  the  w-dependence) 

/  =  ^  UognTn,n'  COS  n'0p' 
n,n' 

+  4  Zl^l  (^"+1  +  9n-l  +g~n+l  +  g-n-\)f^y  COS  n(f)p  COS  Tl' (j)p> 

n,n' 

~4         ^1  (^"-1  +  9-n+l  -  9n+l  -  g-n-l)  f°n'  SlU  n(l)p  Sm  u' (j)p-  .  (13.16) 

n,n' 

This  expression  allows  me  to  write  a  general  system  of  (inhomogeneous)  linear 
equations  for  the  f,''J°,.  In  the  spirit  of  the  choice  for  the  Ui  I  will  assume  that  all 
^,n°  are  negligible  for  n,n'  >  1.  Such  a  choice  is  always  a  solution,  though  it  may 
not  be  the  only  one  in  a  general  case.  Then  only  Qo  and  g2  enter  the  equations 
which  read 

fl,  =  U,  +  U,{g,  +  g,)fl, 
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(13.17) 


where  I  used  the  T!,^  „,  =  T„%, ,  T^„  „,  =  -T„%,.  All  equations  for  indices 
n  ^  n'  have  vanishing  solutions  since  there  is  no  inhomogeneous  term.  Sur- 
prisingly, these  equations  show  a  decoupling  of  partial  waves,  s-wave  scattering 
only  influences  the  s-wave  coefl5cient  of  the  T-matrix  and  similarly  for  rf-wave 
coefl!icients. 

The  solutions  to  the  equations  13.17  are  readily  obtained.  For  example,  I  write 
=  K'^o  +  tlTi  (the  other  r's  are  not  involved).  This  leads  to  two  equations 
for  tl  and  tl  which  I  write  as 


(13.18) 


+  5d 

92  9o 

(A 

^9o  92  ^ 

After  a  bit  of  linear  algebra  I  find 


(A 

[ 

1 

det 

\ 

1  -  ^d92  5d9o 
Sd9o    1  -  Sdg2  j 

with  det  =  (1  -  5,102)^  -  {5d9of-  Finally,  this  leads  to 


(13.19) 


fe    ^  ^dgpT-Q  +  <^d(l  -  5d92)ri 

(1  -  Sd92y  -  {Sd9or 


(13.20) 


Similarly,  I  find 


(13.21) 
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At  q  —  0,  (l)p  =  (f>pi  =  (f) I.  and  the  T-matrix  becomes 

t^,  (oj)  =  foiu)  +  f-{u)  cos  2(1)1,  +  f+{uj)  (13.22) 


with 


{Uo  +  S'i)goro  +  U0T3  +  6sTi 


1  -  (U'o  +  Wo 

,(1- 

,(1- 


(13.23) 
(13.24) 
(13.25) 


The  T-matrix  has  now  an  explicitly  A;-dependent  and  off-diagonal  term  .  Also 
observe,  that  there  is  actually  a  contribution  {f^{u)))  to  the  ifc-independent, 
diagonal  part  involving  the  new  scattering  parameter  5^.  Both  new  contributions 
have  an  impact  on  the  impurity  averaged  self  energy  E,  (see  next  chapter). 

13.3    Constitutive  Equation  for  the  Order  Parameter  Scattering  Strength 

For  simplicity,  let  me  consider  Ss  =  0  first,  that  is,  I  consider  a  d-wave  SC 
without  any  (not  even  a  repulsive)  s-wave  component  in  the  pairing  potential 
V.  This  is  probably  unrealistic,  however,  recent  numerical  work  [89]  has  shown 
that  the  induced  (extended)  s-wave  order  parameter  has  a  vanishing  9  =  0- 
component  due  to  sign  changes  on  the  lattice  diagonals,  and  is  thus  not  relevant 
for  this  discussion. 
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As  before,  I  insert  the  new  T-matrix  into  Eq.  12.5  and  separate  terms.  I  find 


NoV 


(1  -  S,g2) 


Uo9o       ,  Sjgo 


(2008^(20)  J 


(a;2  +  e  +  A2)2 


+ 


where  {..  .)<p  stands  for  the  angular  integral  /  •  •  ••  After  integration  over  the 
energy  ^  I  eUminate  the  l/(A^oV')-term  by  use  of  the  gap  equation  10.6  and  obtain 
(after  canceling  an  overall  constant  27r) 


T  Eu,  (cos'(2(^) 


1  -  Sdg2 


-  S,g2r  -  {S^goY  (a;2  +  A2)3/2 
Sdg2  -  iSdg2f  +  {Sdgo)^        1  \ 


+ 


1  -  [UogoY     (1  -  W  -  +  A2) 


This  lengthy  expression  reduces  properly  to  Eq.  12.11  if  I  expand  the  denom- 
inators and  consider  terms  up  to  first  order  in  5^.  However,  due  to  5^  in  the 
denominators  Eq.  13.27  will  prevent  5^  from  diverging.  Because  of  its  nonlinear- 
ity,  the  solutions  of  this  constitutive  equation  are  only  numerically  accessible.  I 
will  show  the  results  in  the  next  chapter. 

I  am  ending  this  chapter  by  again  pointing  out  the  approximations  used  to 
obtain  the  new  T-niatrix  and  the  constitutive  equation  13.27. 

•  I  assume  short  ranged  order  parameter  perturbations  that  are  characterized 
by  their  volume  integral  5l^k{q  =  0). 
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•  I  consider  a  g-independent  form  of  the  order  parameter  perturbations  with 
a  A;-dependence  motivated  by  the  bulk  order  parameter. 

•  All  the  calculations  are  based  on  the  weak-coupling  BCS-theory. 

I  then  choose  a  pairing  potential  that  limits  the  number  of  partial  waves  con- 
tributing to  scattering.  This  is  not  a  new  approximation,  rather  a  specification 
to  the  case  of  interest.  After  obtaining  the  T-matrix  I  then  use  the  gap  equation 
to  write  a  constitutive  equation  for  the  special  case  of  a  pure  d-wave  pairing 
potential  V.  The  constitutive  equation  allows  the  determination  of  the  only  pa- 
rameter left  in  the  theory,  the  magnitude  of  order  parameter  perturbations  5d. 
Therefore,  is  not  a  free  parameter  but  is  determined  by  the  requirement  of 
self-consistency  of  the  BCS-theory. 


CHAPTER  14 

ORDER  PARAMETER  PERTURBATIONS  AND  SELF  ENERGIES 

14.1    The  Order  Parameter  Scattering  Strength 

The  constitutive  equation  13.27  allows  me  to  self-consistently  determine  the 
magnitude  of  the  order  parameter  perturbations  5d.  In  contrast  to  the  first  order 
result  Eq.  12.11,  from  Eq.  13.27  is  not  just  a  function  of  T/T^  and  Uo  but 
depends  explicitly  on  the  cut-off  frequency  LOg  and  the  coupling  V.  I  choose  for 
the  normal  density  of  states  No  =  0.01,  cOo  =  30  and  V  =  28.3108,  so  that  =  1. 
Therefore,  all  energies  (frequencies,  temperatures)  are  in  units  of  Tc  =  A;b  =  1). 
With  NoV  =  0.283108  I  am  at  the  limit  of  what  is  considered  "weak-coupling". 
Since  A^^  ~  1/Ep,  Ef/T^  ~  100,  appropriate  for  the  High-T^  SC's.  Recall  that 
the  maximum  of  the  bulk  order  parameter  at  T  =  0  is  given  by  A/Tc  =  2.14  for 
a  2D  dj;2_j,2  SC  with  cylindrical  Fermi  surface. 

I  first  show  the  dependence  of  6^  as  a  function  of  the  impurity  scattering 
strength  Uo  at  a  low  temperature  T/T^  =  0.01,  see  Fig.  14.1.  Both  the  first  order 
and  the  fully  self-consistent  (involving  all  orders)  results  rise  monotonically  and 
saturate  as  one  approaches  the  unitary  limit  No-kUo  >  1.  A  peculiar  feature 
is  that  the  all  order  result  is  actually  slightly  higher  than  the  first  order  result 
in  the  Born-limit  NonUo  <  1.  For  large  C/„,  however,  the  first  order  result  is 
higher.  This  I  expected  from  the  fact  that,  for  S^NoTi  ~  1,  the  denominators  in 
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Figure  14.1:  Order  parameter  scattering  strength  vs.  normal  scattering  strength. 
At  r  =  O.OlTc  both  the  first  order  and  all  order  result  for  -5d  rise  monotonically 
and  saturate  in  the  unitary  scattering  limit,  No'nUo  »  1.  The  saturated  value  for 
the  all  order  result  is  about  -SdNgTT  ~  1.2,  with  No  =  0.01  in  units  of  inverse  Tc. 

Eq.  13.27  reduce  the  magnitude  of  5d  compared  to  the  first  order  result.  The  all 
order  result  saturates  at  about  -SdNgTr  ~  1.2.  The  zeroth  order  result  (from  Eq. 
12.5  with  the  standard  T-matrix)  is  about  an  order  of  magnitude  smaller  at  this 
temperatures. 

Next,  I  look  at  the  temperature  dependence  of  the  order  parameter  scattering 
strength.  Fig.  14.2.  For  low  temperatures,  there  is  only  a  weak  T-dependence. 
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Figure  14.2:  Order  parameter  scattering  strength  vs.  temperature.  For  low 
temperature  there  is  only  a  weak  T-dependence.  Close  to  Tc,  the  first  order  result 
(dashed  line)  diverges  like  (1— T/Tc)"^/^  whereas  the  all  order  result  vanishes  with 
the  bulk  order  parameter  A  like  (1  -T/Tc)^/^.  Observe  that  the  maximum  of  the 
all  order  curves  moves  to  lower  temperatures  as  Uo  increases.  For  Uo  ->  oo,  the 
maximum  approaches  T  —  0,  and  the  temperature  dependence  of  -S^  resembles 
that  of  the  bulk  order  parameter  A. 

The  all  order  data  show  a  maximum  that  moves  to  lower  temperatures  as  Ug 
increases.  For  Uo  go,  the  maximum  value  is  probably  (within  numerical  uncer- 
tainty) located  at  T  =  0.  Therefore,  the  temperature  dependence  of  -5^  strongly 
resembles  that  of  the  bulk  order  parameter  A.  As  advertised,  the  behavior  of 
the  first  order  result  is  strikingly  different,  especially  at  temperatures  close  to  Tc 
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where  it  diverges  like  (1  -T/Tc)~^/^.  The  all  order  data  vanish  at  Tc  with  the  bulk 
order  parameter  A  like  (1  —  T/TcYl"^ .  Although  my  local  approximation  for  the 
order  parameter  perturbations  can  be  only  justified  for  temperatures  for  which 
the  Ginzburg-Landau  coherence  length  =  Cr=o/(l  -  T/Tc)^''^  is  less  than  the 
average  inter-impurity  distance,  for  the  small  ^t=q  of  the  high-Tc  SC's  and  low 
impurity  concentrations,  the  temperature  at  which  the  approximation  fails  can 
be  as  high  as  0.99  Tc.  Therefore,  my  approximation  is  justified  for  the  systems  of 
interest  at  all  temperatures  but  the  very  small  GL-regime. 


Having  obtained  the  order  parameter  scattering  strength  5^,  I  now  make  use 
of  the  SCTMA  (see  chapter  10)  to  compute  the  disorder  averaged  self-energy  E. 
Eq.  10.15  is  a  self- consistency  equation  for  the  diagonal  (oc  Tg)  and  off-diagonal 
(oc  Ti)  components  of  E,  since  the  Qo  and  g2  of  the  T-matrix  13.22  are  depending 
on  the  self-energy  via  the  disorder-averaged  Green  function 


The  self-consistency  equation  10.15  must  be  solved  numerically.  I  show  the  re- 
sults below. 

For  better  comparison  with  standard  notation  (e.g.  Ref.  [59]),  I  rescale  the 
scattering  strengths  and  the  functions  Qo  and  §2-  Introducing 


14.2    T-matrix  and  the  Disorder  Averaged  Self  Energy 


G{k,u:)  =  {iu  -  era  -  A^n  -  E^) 


-1 


(14.1) 


1 


1 


c  = 


(14.2) 
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Go  —  go/{No7r)  and  G2  —  g2/{NoTc)  and  F  =  ni/{Non)  (the  scattering  rate  in  the 
normal  state,  i.  e.,  for  T  >  Tc)  the  self-consistency  equation  for  the  considered 
case  reads 

For  a  particle-hole  symmetric  system,  the  rs-component  of  the  self-energy  is  van- 
ishing [93],  so  only  To-  and  ri -components  (S^  and  Ej,  correspondingly)  need  to 
be  determined.  For  real  frequencies  u,  these  self-energy  components  are  complex 
functions  of  In  the  standard  case  without  order  parameter  scattering,  only  Eo 
is  nonvanishing  (as  can  be  seen  from  above  by  setting  =  0).  Then  one  can 
relate  the  imaginary  part  of  the  self  energy  to  a  scattering  rate  or  the  relaxation 
time  (e.g.  see  [59]).  In  the  case  here,  however,  I  also  have  a  ri-component  Ei  of 
the  self  energy  which  also  has  an  imaginary  part.  The  concept  of  a  scattering  rate 
(on  this  microscopic  level)  becomes  now  questionable.  What  should  one  think  of 
a  scattering  rate  that  is  a  matrix  in  particle-hole  space?  Unfortunately,  I  have 
to  abandon  this  intuitive  concept  for  now,  until  other  methods  (like  a  Boltzmann 
equation  approach)  might  show  what  the  equivalent  of  the  standard  scattering 
time  is  for  a  system  with  order  parameter  scattering. 

In  terms  of  the  dimensionless  parameter  c  one  distinguishes  three  regimes  of 
scattering  strengths,  c  >>  1  corresponds  to  the  Born  limit  of  weak  scattering, 
c  ->  0  is  the  strong  scattering  regime,  and  c  ~  1  is  the  intermediate  scattering 
regime.  It  is  instructive  to  replot  Fig.  14.1  as  a  graph  for  c/  vs.  c,  see  Fig.  14.3. 
The  plot  reveals  several  features: 
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Figure  14.3:  Dimeiisionless  order  parameter  scattering  strength  vs.  dimensionless 
normal  scattering  strength.  In  the  Born  limit,  c  »  1,  -c/  >  c,  so  that  order 
parameter  scattering  is  probably  unimportant.  However,  in  the  intermediate 
scattering  regime,  c  ~  1,  -c/  ~  c,  so  both  scattering  processes  are  important.  In 
the  strong  scattering  regime,  c  0,  -c/  remains  finite,  but  can  be  close  to  2/n 
(see  inset),  which  allows  the  new  self  energy  contributions  to  be  large,  see  text 
for  more  detailed  discussion. 

•  In  the  Born  limit  (c  >  1),  -c/  >  c.  Therefore,  order  parameter  scattering 
is  probably  unimportant,  since  the  denominators  of  the  new  contributions 
to  the  self  energy  Eq.  14.3  are  much  larger  than  that  of  the  standard  term. 

•  In  the  intermediate  scattering  regime  (c  ~  1),  -c/  ~  c,  and  by  the  same 
reasoning  as  above  all  contributions  to  the  self  energy  are  important. 


•  In  the  strong  scattering  regime  (c  ->  0),  -c/  remains  finite.  However,  be- 
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cause  of  the  different  structure  of  the  denominators  of  the  new  contributions 
in  Eq.  14.3,  this  does  not  immediately  render  the  new  terms  unimportant 
compared  to  the  standard  term.  For  a  pure  d-wave  SC  G2{uj  =  0)  =  -2/7r, 
so  that  even  for  -cj  ~  1/1.2,  (c/  -  ^2(0;  =  0))  ~  0.2.  This  illustrates,  that 
the  denominator  structure  enhances  the  effect  of  the  order  parameter  scat- 
tering, though  its  bare  strength  is  less  than  that  of  the  standard  scattering 
term. 

I  am  primarily  interested  in  the  last  (strong  scattering)  regime,  since  the  standard 
approach  has  been  successful  in  explaining  the  power  laws  of  thermodynamic 
quantities  like  the  London  penetration  depth  observed  in  the  high-Tc  SC's.  The 
above  argument  stresses  the  necessity  of  a  numerical,  quantitative  evaluation  of 
the  self  energy  components,  since  there  are  too  many  energy  scales  in  the  system 
to  isolate  a  frequency  regime  in  which  power  law  expansions  of  the  self  energy 
are  valid. 

14.3    Diagonal  Self  Energy  Component  Ep  and  the  Density  of  States 

For  unitary  scattering,  i.e.  c  =  0,  in  Fig.  14.4  I  show  the  imaginary  part  of 
the  diagonal  self  energy      (E  =  EoTo  +  EjTi  cos2(^).  For  illustration  purposes, 
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Figure  14.4:  Imaginary  part  of  the  diagonal  self  energy  component  Eq  vs.  fre- 
quency. For  unitary  scattering  (c  =  0)  the  standard  term  (dashed  line)  shows 
a  peak  at  zero  frequency  and  a  characteristic  1  /w-decay  above  an  energy  scale 
7  ~  VTA.  The  new  term  has  a  wide  peak  at  intermediate  frequencies,  but  is 
small  at  lower  frequencies.  As  a  result,  the  DOS  and  the  resulting  thermodynamic 
properties  have  minor  changes  at  intermediate  temperatures,  T  ~  0.57^. 

I  choose  a  rather  large  normal  state  scattering  rate  F  =  0.05.  The  standard  term 
(dashed  line)  has  a  zero  frequency  peak  and  a  characteristic  l/a;-decay  above 
7  ~  VTA.  7  is  roughly  the  energy  scale  below  which  the  self-consistency  aspect 
of  Eq.  10.15  becomes  important  in  the  standard  theory  without  order  parameter 
scattering.  The  new  term  has  a  wide  peak  at  intermediate  frequencies,  but  is 
small  at  lower  frequencies.  -JmE^  is  the  main  ingredient  in  the  determination 
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Figure  14.5:  Quasiparticle  density  of  states  (DOS)  vs.  frequency.  The  pure  d- 
wave  SC  shows  linear  dependence  at  low  frequencies,  the  standard  theory  (unitary 
scattering,  c  =  0)  of  a  dirty  d-wave  SC  gives  a  basically  flat  DOS  for  frequencies 
below  7  ~  y/TA.  This  translates  in  linear  (pure  d-wave)  and  quadratic  (dirty  d- 
wave)  behavior  for  the  London  penetration  depth,  in  agreement  with  experiments 
on  clean  and  dirty  samples.  The  new  scattering  due  to  Sa  leads  only  to  a  mild 
hump  at  intermediate  frequencies,  even  for  the  rather  large  value  of  F  =  0.05 
(corresponding  to  a  very  impure  sample).  Therefore,  the  new  scattering  has  no 
qualitative  bearing  on  thermodynamic  properties. 

of  the  quasiparticle  DOS,  see  Eqs.  10.7  and  14.1.  Since  the  low  energy  behavior 
of  — /mEo  is  essentially  unchanged,  so  are  the  low  energy  DOS  and  the  resulting 
thermodynamic  quantities.  I  show  the  new  DOS  together  with  the  standard 
result  and  the  DOS  of  a  pure  d-wave  SC  in  Fig.  14.5.  The  new  term  affects 
the  DOS  only  at  intermediate  frequencies  and  the  thermodynamics  only  mildly 
at  intermediate  temperatures,  T  ~  0.5Tc.  Not  just  for  the  sake  of  completeness, 
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Figure  14.6:  Real  part  of  the  diagonal  self  energy  component  vs.  frequency. 
ReEo  is  proportional  to  u  at  frequencies  below  7  ~  vTA.  The  new  term  again 
adds  a  hump  at  intermediate  frequencies.  ReEg  has  usually  only  minor  impact 
on  physical  observables. 

I  also  show  the  real  part  of  in  Fig.  14.6.  ReHg  has  usually  small  impact  on 
physical  quantities.  It  is  important,  though,  that  at  frequencies  below  7  it  rises 
linearly  from  zero,  as  it  does  without  order  parameter  scattering.  Again,  the  new 
term  adds  a  hump  at  intermediate  frequencies. 


14.4    Results  for  the  Off-Diagonal  Self  Energy  Component  Ei 


The  off-diagonal  self  energy  component  Ei  is  an  entirely  new  term,  so  I  can  not 
compare  to  the  standard  theory  here.  Instead,  in  Fig.  14.7  I  show  its  frequency 
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Figure  14.7:  Imaginary  part  of  the  off-diagonal  self  energy  component  Ei  vs. 

frequency.  /mEi  is  finite  and  positive  below  an  c  (5(i)-dependent  energy  scale. 

With  increasing  c  (decreasing  -6d)  ImEi  becomes  smaller  in  overall  magnitude. 

dependence  for  various  scattering  strengths  c.  The  most  important  aspect  of  the 
curves  is  the  fact  that  /m£i  is  finite  and  positive  below  a  new  energy  scale  that 
in  the  strong  scattering  regime  depends  only  weakly  on  c  (or  more  accurate,  on 
6d).  Observe,  that  the  crossover  to  positive  values  of  /mEi  stays  almost  the  same 
for  c  <  1,  and  only  slowly  moves  to  larger  frequencies  as  c  increases  (and  -6d 
decreases).  The  crossover  frequency  does  not  depend  on  7  in  any  simple  fashion, 
and  for  the  reasonable  values  of  the  normal  scattering  rate  the  crossover  frequency 
is  larger  than  7.  On  general  symmetry  grounds  ImEi  must  be  antisymmetric 
about  the  Fermi  level.   The  surprise  is  that  rather  than  going  through  zero, 
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7mEi  has  a  finite  step  at  the  Fermi  level.  This  feature,  in  conjunction  with  the 
linear  behavior  of  Re'Eo,  makes  the  new  self  energy  component  important  when 
it  comes  to  transport  properties,  as  will  be  shown  in  the  next  chapter. 
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Figure  14.8:  Real  part  of  the  off-diagonal  self  energy  component  Ei  vs.  frequency. 
ReT,i  is  finite  and  negative  at  all  frequencies  below  the  bulk  order  parameter  A.  It 
is  a  frequency  dependent  renormalization  of  the  bulk  order  parameter.  However, 
since  it  is  always  much  smaller  than  A  itself,  it  is  of  negligible  importance. 

The  real  part  of  the  off  diagonal  self  energy  component  Ei  is  of  minor  impor- 
tance when  it  comes  to  physical  observables.  It  is  instructive  to  realize,  that  it 
is  ReEi,  rather  than  Sd,  that  renormalizes  the  bulk  order  parameter  A.  Re^i  is 
negative  for  all  frequencies,  and  therefore  it  is  a  frequency  dependent  reduction 
of  the  bulk  order  parameter.  However,  it  is  much  smaller  in  magnitude  than  the 
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bulk  order  parameter  A^.  Since  it  enters  equations  always  in  addition  to  A  it  has 
negligible  effects.  This  justifies  partially  our  use  of  BCS  rather  than  Eliashberg 
theory,  since  large  frequency  dependent  corrections  to  A^  would  necessitate  the 
use  of  the  strong  coupling  Eliashberg  theory. 


CHAPTER  15 
APPLICATION:  MICROWAVE  CONDUCTIVITY 


As  an  example  of  the  effect  of  the  order  parameter  scattering  on  transport 
properties  I  will  discuss  the  "microwave"  conductivity  of  a  d-wave  SC.  It  is  noth- 
ing but  the  conductivity  in  the  case  of  a  extreme  Type  II  SC  (like  the  high-Tc 
SC's)  where  the  penetration  depth  is  much  larger  than  the  coherence  length. 
It  then  suffices  to  consider  wave  vectors  q  of  the  current  response  that  are  ap- 
proximately zero.  The  microwave  conductivity  a{T,  Q)  {fl  being  the  microwave 
frequency)  is  in  general  a  tensor  that  within  linear  response  is  determined  by  the 
retarded  current-current  correlation  function  (Kubo  formula).  A  detailed  discus- 
sion of  the  formalism  can  be  found  in  Refs.  [59,94].  I  just  quote  the  result  for 
the  diagonal  components  of  the  conductivity  tensor. 

a{T,  n)  =  ——  /     duj  [tanh(/Ja;/2)  -  tanh(/?(a;  -  fi)/2)]  5(a;,  Q,  T),  (15.1) 

Z'III\l  J —  oc 

where  n  is  the  electron  density,  e  and  m  the  electron  charge  and  mass,  respectively, 
and 


5(a;,Q,r)  =  -ilm/-#  ( ^'4^^  +  ^'^)  +  ^4^,^  -  ^',^) 
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Here,  Cj±  —  (^[ixi  rk  zO"*")  with  w  =  Ci;  —  So(a;),  A^^^  =  A^(A,^  ±  zO"*")  with  A^  = 
A^4-Si(a;)  cos  2^,  and  ^±  =  %<gTi[ijj)\jib\  —  A^  _,_.  The  infinitesimal  izO"*"  indicates 
whether  one  has  to  evaluate  the  self  energies  above  or  below  the  branch  cut  in 
the  complex  w-plane.  The  primes  indicate  that  the  corresponding  quantities  are 
to  be  evaluated  at  w  —  ^7. 

In  the  limit  of  vanishing  microwave  frequency  Q  — >  0  the  expressions  can 
be  somewhat  simplified.  For  T  >  f2,  [tanh(/?a;/2)  -  tanh(/?(a;  -  Q)/2)]/2ri  -)■ 
—dY{ijj)lduj  with  ¥{yj)  the  Fermi  function.  The  first  term  in  the  integrand  of 
S{u},Q,,T)  becomes  simply  Co\/(^.  This  term,  however,  turns  out  to  give  only 
a  very  small  contribution  at  any  experimentally  accessible  temperature.  Nev- 
ertheless, it  has  theoretical  significance,  since  at  T  =  0  it  gives  a  nonzero  con- 
tribution to  the  conductivity  that  is  independent  of  the  impurity  concentration, 
CTo  =  ne^/(m7rA),  as  first  pointed  out  by  Lee  [95].  The  second  term  in  the  inte- 
grand of  Eq.  15.2  remains  complicated,  but  it  is  instructive  to  write  it  partly  in 
terms  of  the  real  and  imaginary  parts  of  the  self  energies.  Multiplying  numerator 
and  denominator  with  (^+  +       I  obtain  (ignoring  an  overall  factor) 

{u  -  ReT,o)u-  -  /mSi  cos20A^+  Re(+ 
-ImEoiu  -  ReEo)  -  /mEi(A  +  i?eEi)cos2(2(/))  ^^^"^^ 

For  small  frequencies  I  showed  in  the  previous  chapter  that  ImEo  is  finite  and 
negative,  ReEo  vanishes  linearly  with  lo,  however  /mEi  is  nonzero  and  positive. 
This  implies  that  the  first  fraction  of  15.3  has  a  pole  at  a  certain  angle  (f)^.  This 
pole,  however,  is  countered  by  a  zero  in  Re^+.  Together,  the  term  15.3  has  a  cusp 
at  the  angle  (po  for  small  frequencies.  Such  features  do  not  exist  in  the  standard 
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theory  with  no  Ej  present.  Thus,  this  analysis  hints  to  the  importance  of  the 
new  off-diagonal  self  energy  Sj. 


Figure  15.1:  Microwave  conductivity  vs.  temperature.  For  the  normal  state 
scattering  strength  F  =  0.001  I  show  the  cr(T,  Q  =  0)  for  the  new  theory  including 
order  parameter  scattering  and  the  standard  theory,  both  in  the  unitary  scattering 
limit  c  =  0.  At  very  low  temperatures  the  curves  are  almost  identical.  Above 
T  ~  O.lTc  the  new  theory  gives  a  much  lower  conductivity  due  to  the  hump  in 
Im^o-  The  graph  is  roughly  linear  until  it  reaches  a  maximum  just  below  Tg. 
At  Tc  the  conductivity  is  given  by  a{T^)  =  The  standard  theory  has  a 

much  broader  'coherence'  peak  until  it  also  reaches  the  normal  state  conductivity 
a(Te). 

In  Fig.  15.1  1  show  the  zero-frequency  microwave  conductivity  a(T,  =  0) 
for  the  new  theory  including  order  parameter  scattering  compared  to  the  stan- 
dard theory  without  order  parameter  scattering.  At  very  low  temperatures  the 
curves  are  almost  identical  and  behave  hke  T^.  Above  T  ~  O.lTc  the  new  the- 
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ory  has  much  lower  conductivity  due  to  the  hump  in  ImHo  which  leads  to  an 
increased  scattering  rate  compared  to  the  standard  theory  without  order  param- 
eter scattering.  The  graph  of  the  new  theory  is  roughly  linear  until  it  reaches  a 
maximum  just  below  Tc.  At  Tc  the  conductivity  should  be  given  by  the  normal 
state  conductivity,  (j{Tc)  =  The  new  theory  fulfills  this  requirement  due 

to  the  vanishing  of  5^  at  Tc.  The  standard  theory  has  a  much  broader  peak  (due 
to  the  'coherence  factors'  in  the  language  of  BCS-theory).  It  rises  much  faster  at 
intermediate  temperatures  and  reaches  the  normal  state  conductivity  o'(Tc),  too. 

In  Fig.  15.2  I  show  the  low  temperature  part  of  Fig.  15.1.  As  mentioned 
above,  below  T  ^  O.lTc  there  are  only  small  differences  between  the  new  theory 
including  order  parameter  scattering  and  the  standard  theory.  The  new  theory 
has  a  slightly  enhanced  conductivity  above  T  ~  7  ~  y/TA  which  can  be  traced  to 
the  new  self  energy  contribution  Ei.  This  enhancement  leads  to  approximately 
linear  T-dependence  above  T  ~  7.  However,  below  that  temperature  the  conduc- 
tivity behaves  like  as  it  did  in  the  standard  theory  without  order  parameter 
scattering.  Above  T  ~  O.lTc  the  hump  in  /mE^  leads  to  a  suppression  of  the 
conductivity  due  to  the  increased  scattering.  The  overall  shape  of  the  conductiv- 
ity of  the  new  theory  appears  linear  after  the  initial  quadratic  rise.  In  contrast, 
the  standard  theory  is  essentially  quadratic  over  the  total  temperature  range  of 
the  graph. 

I  mentioned  in  the  previous  chapter  that  the  order  parameter  scattering  might 
also  be  of  importance  in  the  case  of  intermediate  normal  scattering  strength. 
For  c  =  0.3  (which  corresponds  to  Uo  of  about  the  Fermi  energy)  I  show  again 
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Figure  15.2:  Microwave  conductivity  vs.  temperature  at  low  temperatures.  For 
the  normal  state  scattering  strength  F  =  0.001  and  in  the  unitary  scattering 
limit  c  =  0  the  new  theory  including  order  parameter  scattering  and  the  standard 
theory  show  only  small  differences  below  T  ~  O.lTg.  The  new  theory  gives  slightly 
enhanced  conductivity  but  still  behaves  like  at  very  low  temperatures  T  <  7. 
Above  T  ~  O.lTc  the  new  theory  gives  a  much  lower  conductivity  due  to  the 
hump  in  ImLo-  The  graph  of  the  new  theory  appears  linear  for  T  >  7.  The 
standard  theory  behaves  like  T"^  for  the  total  temperature  range  in  this  graph. 

the  zero-frequency  microwave  conductivity  for  both  the  standard  and  the  new 
theory.  As  before  there  is  an  enhancement  at  low  temperatures  and  suppression  of 
the  conductivity  at  higher  temperatures.  The  low-T  enhancement  is  much  more 
visible  than  in  the  unitary  scattering  case.  However,  the  major  difference  between 
the  curves  is  still  at  intermediate  temperatures.  The  large  residual  conductivity 
at  low  T  is  in  strong  contrast  to  experiments  on  clean  high-Tc  samples.  Therefore, 
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this  case  is  of  pure  theoretical  interest  for  the  time  being,  lacking  an  experimental 
realization  at  the  moment. 
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Figure  15.3:  Microwave  conductivity  vs.  temperature  for  nonunitary  normal 
scattering.  For  F  =  0.001  I  show  the  g{T,  =  0)  for  the  new  theory  and  the 
standard  theory,  for  an  intermediate  scattering  strength,  c  =  0.3.  At  low  temper- 
atures the  new  theory  gives  a  slightly  higher  conductivity.  However,  the  major 
deviation  is  at  higher  temperatures,  where  the  conductivity  is  suppressed  due  to 
the  additional  hump  in  IniLo- 


So  far,  I  have  only  considered  the  case  of  vanishing  microwave  frequency, 
1}  =  0.  In  the  standard  theory  finite  Q  lower  the  conductivity,  but  do  not  change 
the  low  temperature  power  law.  I  do  not  expect  this  to  be  different  in  the  theory 
including  order  parameter  scattering.  Inelastic  scattering,  e.g.  due  to  spin  fluc- 
tuations, is  negligible  at  low  temperatures  but  rises  approximately  like  {T/T^f, 
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and  it  therefore  becomes  important  at  intermediate  temperatures  (depending  on 
r).  Beyond  such  a  crossover  temperature  inelastic  scattering  dominates,  hence 
the  conductivity  will  decrease.  The  increase  at  low  T  due  to  impurity  scattering 
and  the  decrease  at  higher  T  due  to  inelastic  scattering  give  rise  to  the  peak  in 
the  conductivity  seen  in  experiments.  This  peak  is  thus  not  related  to  the  'co- 
herence' peak  shown  in  Fig.  15.1.  I  have  performed  preliminary  calculations  on 
the  conductivity  including  finite  frequencies  and  inelastic  scattering.  However, 
because  of  their  preliminary  nature  I  refrain  from  including  these  results  in  this 
thesis. 


CHAPTER  16 
CONCLUSIONS  TO  PART  II 


I  have  investigated  the  effect  of  nonmagnetic  impurities  on  a  d-wave  SC.  In 
addition  to  the  standard  potential  scattering,  I  included  the  scattering  due  to 
local  perturbations  of  the  order  parameter  around  the  impurity  site.  Since  these 
perturbations  are  not  a  priori  small,  I  determined  them  self-consistently  via  the 
gap  equation.  To  do  that,  I  employed  a  local  approximation  in  which  I  replaced 
the  short  ranged  spatial  dependence  of  the  perturbations  by  a  Dirac  5-function. 
I  was  able  to  analytically  compute  the  momentum  dependent  T-matrix  of  a 
single  impurity  in  a  bulk  d-wave  SC.  I  analyzed  the  validity  of  the  standard  self- 
consistent  T-matrix  approximation  (SCTMA)  for  the  case  of  a  2D  d-wave  SC, 
and  concluded  that  the  SCTMA  is  qualitatively  sufficient  to  describe  the  effects 
of  dilute  concentrations  of  impurities  on  such  systems.  I  numerically  determined 
the  disorder  averaged  self  energy  E  of  the  bulk  system  employing  the  SCTMA. 

Two  new  contributions  show  up  as  a  consequence  of  the  new  scattering  mech- 
anism, one  modifying  the  diagonal  self  energy  Eo,  the  other  one  a  completely 
new  off-diagonal  term  Si  with  the  same  angle  dependence  as  the  bulk  order 
parameter.  The  new  self  energy  contributions  have  negligible  effect  on  thermo- 
dynamic properties,  since  they  add  only  a  minor  feature  to  the  density  of  states 
away  from  the  Fermi  level.  Their  impact  on  transport  properties  is  more  signifi- 
cant. At  low  temperatures  the  off-diagonal  contribution  Si  slightly  enhances  the 
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microwave  conductivity  (at  zero  frequency),  making  it  appear  linear  for  tempera- 
tures T  >  7  ~  \/rA.  However,  this  effect  is  small,  and  the  very  low  temperature 
behavior  ~  is  unchanged.  At  higher  temperatures  the  second  contribution 
has  a  more  obvious  impact,  since  it  reduces  the  conductivity  significantly  in  com- 
parison to  the  standard  theory  without  order  parameter  scattering.  However, 
because  of  the  dramatic  increase  of  inelastic  scattering  at  higher  temperatures 
(not  considered  this  far)  this  suppression  of  the  conductivity  will  be  experimen- 
tally obscured  due  to  the  much  stronger  suppression  as  a  consequence  of  inelastic 
scattering. 

It  is  thus  questionable,  whether  experiments  of  the  current  resolution  will  be 
able  to  detect  the  effects  of  order  parameter  scattering  on  transport  properties. 
In  the  real  high-Tg  materials  there  might  also  be  other  scattering  mechanisms 
due  to  other  impurities  than  considered.  From  the  theorists  point  of  view,  the 
new  theory  can  not  explain  the  linear  temperature  dependence  at  very  low  T, 
though  it  provides  a  step  in  the  right  direction.  Open  questions  concern  the 
dependence  of  the  conductivity  on  the  external  frequency  Q  and  the  influence  of 
inelastic  scattering. 

These  issues,  however,  will  be  addressed  at  a  later  time  and  not  within  the 
framework  of  this  thesis.  For  the  time  being  I  will  rest  my  case. 


CHAPTER  17 
FINAL  CONCLUSIONS 


I  considered  t  wo  systems  with  small  amounts  of  impurities  added  to  an  other- 
wise pure  material.  In  tiie  first  part  I  studied  the  conductance  and  susceptibility 
of  clean  metal  point  contacts  and  metal-insulator-metal  tunnel  junctions  in  the 
presence  of  impurities  with  two  states  of  degenerate  energy.  Depending  on  the 
kind  of  interaction  of  the  impurity  with  the  conduction  electrons,  the  low  temper- 
ature behavior  of  the  physical  properties  was  qualitatively  changed,  compared  to 
a  system  with  standard  potential  scattering.  I  have  studied  the  nonequilibrium 
aspects  of  these  systems  in  the  case  of  a  static  bias  V.  Scaling  of  the  nonlinear 
conductance  was  observed  at  low  temperatures  as  well  as  deviations  from  scaling 
at  higher  temperatures. 

In  the  second  part  I  studied  the  influence  of  standard,  nonmagnetic  impu- 
rities on  a  d-wave  superconductor.  The  nodes  of  the  order  parameter  of  the 
d-wave  SC  lead  to  qualitative  different  behavior  for  both  thermodynamical  and 
transport  properties  when  compared  to  the  metallic,  s-wave  SC's.  I  included 
the  response  of  the  order  parameter  to  the  impurity  potential  as  an  additional 
scattering  mechanism.  Though  important  at  intermediate  temperatures,  the  low 
temperature  behavior  of  the  microwave  conductivity  was  unchanged  when  com- 
pared to  standard  theories. 

In  conclusion,  the  first  system  showed  dramatic  qualitative  changes  of  observ- 
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ables  upon  addition  of  small  amounts  of  impurities.  Such  changes  are  present 
also  in  the  second  system,  however,  my  modification  of  the  standard  treatment 
turned  out  to  be  a  quantitative  effect,  leaving  the  qualitative  behavior  unaffected. 
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APPENDIX  A 
ENFORCING  THE  CONSTRAINT 


During  the  iteration  procedure  the  lesser  function  can  change  quite  dramat- 
ically if  the  Ijias  and/or  temperature  is  changed  substantially.  It  is  not  obvious 
that  these  changes  lead  to  functions  that  still  fulfill  the  constraint 

Z  =  I  de[Na{e)  +  Mb{e)]  =  1  .  (A.l) 

I  would  like  to  enforce  the  constraint  without  disturbing  the  result  of  the  iteration. 
It  does  not  matter  whether  I  modify  the  algorithm  during  the  iterations,  as  long 
as  the  modifications  vanish  when  convergence  is  achieved.  I  introduce  a  constant 
Xo  in  a  way  that  allows  me  to  satisfy  the  constraint  after  each  iteration.  Defining 

^         _   E<(a;)  

"'^"^'^^  -  (uj-e,  +  X,-  ReE'-(w))2  +  (ImE'-(a;))2  ^^'^^ 

'^"^^f  ~  ^uj  +  X,-  Ren^{uj)y  +  (Imn'-(a;))2  ^^'^^ 
I  determine  A„  by  the  requirement  of  fulfilling  the  constraint 

I  d€[Nax^  (e)  +  Mb^,^  {e)]  =  1  .  (A.4) 

A  solution  can  be  found  by  standard  methods  since  the  integral  is  a  monotonically 
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decreasing  function  of  Aq.  After  finding  Ao  I  define  new  spectral  functions  via 

'  (a;-e^  +  Ao-ReE'-(a;))2  +  (Im£'-H)2  ^'^^ 

,  -  _  Imn''(a;) 
'"^^^  ~  [oj  +  Xo-  Ren'-(a;))2  +  (Imn'-(a;))2  '  ^ 

I  then  use  these  four  functions  as  input  for  the  next  iteration.  Observe  that  if  I 
introduced  Ao  also  in  the  arguments  of  the  self  energies,  I  would  just  have  shifted 
the  energy  scale  by  A^.  In  that  case  Ao  would  be  a  modification  of  the  chemical 
potential  ft.  During  the  iterations,  when  the  functions  and  Ao  are  constantly 
changing,  the  introduction  of  Ao  is  a  significant  modification  of  the  'true'  spectral 
functions.  However,  as  one  approaches  convergence  the  changes  in  the  lesser  and 
spectral  functions  become  smaller  and  smaller,  and  so  are  the  changes  in  Ao.  At 
the  point  of  convergence  the  changes  are  negligible  and  I  arrive  at  functions  that 
are  solutions  to  the;  NCA-equations  and  fulfill  the  constraint.  The  iteration  pro- 
cedure incorporated  A^  into  the  solutions,  so  that  at  convergence  Ao  is  indeed  just 
a  shift  in  the  encrg}'  scale,  though  it  was  originally  introduced  rather  artificially 
to  enforce  the  constraint. 

It  turns  out  that  the  shift  due  to  Ao  has  another  effect  that  is  quite  important 
for  the  numerical  procedure.  As  will  be  explained  in  more  detail  in  Appendix  B 
it  is  crucial  to  use  meshes  that  provide  high  resolution  of  all  sharp  features  of  the 
integrand.  The  lesser  and  spectral  functions  have  sharp  peaks  which  would  shift 
during  the  iteration  procedure.  Thus,  one  would  have  to  modify  the  mesh  and 
interpolate  the  functions  after  each  iteration.  However,  the  changing  Ao  com- 
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pensates  these  shifts  and  allows  me  to  work  with  a  fixed  mesh.  This  saves  a 
considerable  amount  of  CPU-time,  which  is  partially  the  reason  why  I  can  reach 
much  lower  temperatures  than  previously  obtained. 


APPENDIX  B 

INTEGRATION  MESHES  FOR  EQUILIBRIUM  AND  NONEQUILIBRIUM 

NCA 

In  principle,  tlie  integrals  in  the  NCA-equations  run  from  — oo  to  oo  since  I 
am  using  a  gaussian  density  of  states.  The  smallest  features  of  the  integrands 
have  a  width  given  by  the  temperature  T  which  can  be  of  the  order  lO"*^  in  terms 
of  the  width  of  the  gaussian.  Therefore,  a  linear  (equidistant)  integration  mesh 
with  spacing  of  about  T  is  out  of  the  question;  it  would  be  too  time  consuming. 
However,  there  are  only  few  features  of  that  width  whereas  most  of  the  integrand 
has  structure  on  a  scale  determined  by  bare  level  width  F.  So  I  need  to  construct 
a  mesh  that  achieves  the  resolution  of  all  small  scale  features  like  the  spectral 
function  peaks  and  tlie  Fermi  (distribution)  function  step(s)  without  wasting  lots 
of  mesh  points  in  regions  where  the  resolution  can  be  much  more  coarse. 

In  equilibrium,  tliis  can  be  achieved  by  defining  the  mesh  spacing  far  from  the 
Fermi  level  (ep  =  0)  according  to  a  inverse  tangent-function  up  to  an  interface 
point  w/.  There  it  is  matched  with  a  logarithmic  mesh  that  leads  up  to  the  Fermi 
level.  In  principle,  this  is  nothing  but  a  transformation  of  the  integration  variable 
which  I  choose  corresponding  to  my  needs  in  different  regions  of  the  frequency- 
axis.  Consider  L  equidistant  meshes  {x{},  i  =  1. .  .m  oi  m,  I  =  1 . . .  L  points.  I 
map  these  meshes  to  nonuniform  meshes  {a;-}  via  functions  h'{x')  : 

uj[^h'{x\),     i  =  l...n'  (B.l) 
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A.;  =  ^Ax'  (B.2) 


Here,  Ax'  are  ilic  .si)a(  iiigs  of  the  equidistant  meshes.  I  can  now  rewrite  the 
numerical  integration  of  an  arbitrary  function  k{u!)  over  the  nonuniform  meshes 
{ujI}  as  an  integration  over  the  equidistant  meshes  {x\}: 


2^1     —  l^i^i 


+  ^MO)  .  (B.3) 

The  a}  =  uj[,b^  =  uj\^^  are  the  limits  of  integration  of  the  different  regions  of  the 
frequency-axis  .  To  cover  the  whole  axis  I  must  have  a'"*"^  =6'.  In  equilibrium, 
I  can  get  by  with  four  regions:  [-oo,  -w/),  =  e^),  [0,a;/),  [a;/,oo],  where 

u)i  is  an  interface  frequency  (le^^l  -  T  >  uj  »  T^)-  If  I  choose  the  functions 
h}{x'')  as  tan(x')  in  the  regions  with  large  absolute  frequency  and  as  exp(a;')  in 
the  regions  |a;|  <  uj[  \  can  create  large  mesh  point  spacings  far  from  tp  and 
exponentially  small  spacings  ('logarithmic'  mesh)  at  tp  =  0.  Proper  adjustments 
of  constants  in  the  /)/'s  is  required  to  ensure  a  smooth  crossover  of  the  mesh  point 
spacing  at  the  interface  points  and  a  minimal  spacing  at  e^-  of  at  least  10  times 
smaller  than  T  (and/or  V  out  of  equilibrium). 

Crucial  for  the  success  of  this  procedure  is  the  introduction  of  (see  Ap- 
pendix A)  in  the  iteration  procedure.      shifts  the  peaks  of  the  spectral  functions 
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to  the  neighborhood  of  ei?  and  keeps  them  there  at  every  iteration.  Without  Aq 
the  peak  positions  would  shift.  So  even  if  I  started  out  with  good  resolution  of 
the  peaks  the  iteration  procedure  would  move  them  to  regions  with  worse  reso- 
lution. This  would  force  me  to  adjust  the  integration  mesh  after  each  iteration 
and  interpolate  the  functions  on  the  new  mesh  at  great  computational  cost  in 
order  to  keep  the  resolution.  Thus,  the  introduction  of  Aq  kills  two  birds  with 
one  stone:  It  enforces  the  constraint  and  allows  good  resolution  of  all  features 
without  changing  the  integration  mesh  after  every  iteration. 

Out  of  equilibrium  the  distribution  function  is  a  double  step  function  with 
steps  at  ±V/2.  It  tui  ns  out  that  (in  the  Kondo  limit)  the  slave  boson  spectral  and 
lesser  functions  show  (broadened)  peaks  at  about  the  same  frequencies.  However, 
the  pseudofermiou  functions  do  not  behave  as  nicely.  They  do  not  split,  but  have 
a  single  peak  somewhere  between  the  Fermi  level  and  V/2  that  shifts  not  linearly 
with  V.  To  cope  with  such  behavior  I  wish  to  have  good  resolution  at  ±V/2  and 
at  cf-  The  latter  one  is  to  improve  the  resolution  at  the  location  of  the  peak  of 
the  pseudofermiou  functions.  Unfortunately,  I  do  not  know  how  this  location  will 
move  with  increasing  V. 

To  achieve  good  resolution  at  the  mentioned  frequencies  I  let  the  logarithmic 
mesh  end  at  ±V/2  (c  oming  from  larger/smaller  frequencies)  and  choose  the  spac- 
ing in  between  according  to  the  sum  of  two  tanh-functions  which  have  their  zero 
shifted  to  ±V/4,  respectively.  I  have  to  choose  parameters  of  these  functions, 
so  that  the  mesh  spacings  at  the  crucial  energies  is  small  enough  to  resolve  all 
features  of  the  integrand.  These  parameters  depend  on  the  bias  V.  They  have 
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to  be  calculated  befoie  the  mesh  is  'set  up'  whenever  I  change  the  potential  from 
one  run  to  the  next.  However,  once  the  mesh  is  set  I  do  not  have  to  change  it 
anymore  during  the  iteration,  because  of  the  same  reasons  as  in  equilibrium. 

The  minimal  number  of  integration  points  is  200  and  250  for  equilibrium  and 
out  of  equilibrium,  respectively.  Out  of  equilibrium  I  need  about  50  points  more 
for  the  'inner'  region  between  ±V/2  at  moderate  bias  V  <  2QTk-  For  higher 
bias  I  have  to  introduce  more  points  in  the  inner  region.  Convergence  is  achieved 
within  100  -  200  iterations.  The  total  running  time  on  a  current  workstation  is 
less  than  5  minutes  for  equilibrium  and  about  10  minutes  out  of  equilibrium. 


APPENDIX  C 

g-DEPENDENT  ORDER  PARAMETER  PERTURBATIONS  TO  FIRST 

ORDER 

Here,  I  compute-  tlie  g  dependence  of  small  order  parameter  perturbations 
Sd{q)  in  a  d-wave  SC.  I  use  the  gap  equation  with  first  order  self-consistency  with 
respect  to  S^iQ)-  Fiuthermore,  I  assume  that  q  is  small  compared  to  the  Fermi 
momentum  A;. .  Tlicn  all  important  g-dependence  comes  from  the  ^-part  of  the 
bulk  Green  functions  (not  from  A^.)  I  obtain  in  straightforward  generalization  of 
the  procedure  outlined  in  Chapter  12 


1.^ 


SM  il/V  +  r.V„  ^  I  ^Cos\2cj>)  [(AJ  -     +  \{v,q)')IM)  -  h{ct>) 


-TiV^Ey  iCos^(2^)^3|^A..,(,) 


(C.l) 


where  the  integrals  A (^),  72(0)  are  given  by 


IM  =  / 


W)  =  / 


e 

[u^  +  A^  +     +  \v,q)^){u^^  +       +  (e  -  \vpq)^)  ^^'^^ 

Here,  I  also  ignored  all  r/  dependence  which  would  arise  from  a  nonlinear  band 
about  the  Fermi  le\  el.  This  is  justified  for  small  q  and  the  Fermi  energy  not  to 
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close  to  the  band  edge.  The  ^  integrals  can  be  performed  and  give 


TT 


^1^^)  -         ,2   ,    A2M/2,  .2    ,    A2    ,    W,-r    -.2  (^'^^ 


As  for  the  =  0  case  I  can  eliminate  the  1/l^-term  in  Eq.  C.l.  The  final  equation 
for  5d.{q)  reads 


A  similar  result  was  previously  derived  by  Choi  [90]  by  means  of  a  quasi-classical 
approximation.  However,  his  result  has  either  a  misprint  or  is  only  valid  in  the 
Ginzburg-Landau  regime,  i.e.  in  the  case  of  a  small  order  parameter  A^.  The 
above  result  is  valid  for  all  temperatures  under  the  above  mentioned  conditions. 
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